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III — La convection vive d’un fluide actif 
s'écoulant en régime laminaire permanent 
le long d'une plaque (suite). 


6. Interprétation intuitive du coefficient de 
convection laminaire vive. — On voit sur les for- 
mules (33) à (35) que le coefficient de convection 
se scinde en deux termes : 

Le premier, (x:),, qui reste seul quand mQ}; — 0, 
coïncide avec le coefficient classique de Pohlhausen 
fourni par les formules (22) et (23). Nous l’appel- 
lerons coefficient de convection morte, pour marquer 
qu'isolé, il suppose un fluide entièrement passif, 
véhiculant une chaleur qu’il ne contribue pas à 
créer ou à détruire. 

Le second terme, («.),, est nouveau. Nous l’appel- 
lerons terme vif, pour rappeler qu’il exprime l'effet 
supplémentaire de la chaleur dégagée ou détruite 
par le fluide, dans une flamme par exemple. 

On voit tout de suite que le rapport du terme vif 
au coefficient de convection morte de Pohlhausen 
est égal à 


mQ 
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La présence de tous les termes qui figurent dans 
cette expression s'explique intuitivement : 

Si l'écart thermique T, — 0, est grand, le méca- 
nisme classique de la Cnrectibn morte est actif, 
et il est naturel que l’effet du réchauffage interne 
de la flamme soit moins sensible. 

Même remarque si l’on a affaire à un fluide de 
grosse capacité calorifique co, par unité de volume, 
car il véhicule beaucoup de chaleur, et le fait qu’il 
brûle en ajoute relativement moins. On a vu que Co 
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(Eve) 


; : de à FE 
résulte du produit Cp où disparaît l'effet de la 
0 


dilatation sur le nombre des molécules par unité 
de volume, tout au moins pour une réaction équimo- 
léculaire; pour une réaction qui modifierait beaucoup 
le nombre des molécules en cause, la question serait . 
à reconsidérer, mais l'approximation faite est suffi- 
sante pour les flammes. 

On notera que l’annulation de l'écart thermique 
total T. — 1), ou de la capacité calorifique ce, rend 
infini le terme vif du coefficient de convection; nous 
y reviendrons. 

Pour la chaleur totale Q; apportée par la combus- 
tion dans l'unité de volume, il est compréhensible 
que le terme vif croisse dans le même rapport qu’elle. 

Le facteur m, qui est en somme le coefficient 


d'efficacité de Q@,, passe comme on l’a vu par un 


maximum très plat quand la combustion est en 
pleine évolution. Il exprime au fond que le fluide 
réagit effectivement, par exemple que les gaz de la 
flamme brûlent vraiment; on conçoit bien en effet 
que quand la réaction n’est pas encore commencée 
ou qu’elle est terminée, le profil T(y) des tempé- 
ratures à la paroi se redresse aussitôt pour venir se 
confondre avec celui de la convection morte (fig. 3) 

Le rapport entre les coefficients de convection vraie 
(ou totale) et morte est représenté par l’ensemble 
des termes entre crochets des formules (33) à (35) 
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Il mesure la proportion dans laquelle le change- 
ment du profil T(y) des températures dû à la réaction 
affecte le transfert de chaleur, pour un même écart 
thermique total T, — 4. 
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7. Discussion du coefficient de convection 
laminaire vive. 


1° Le signe du terme vif (x:), peut changer : positif 
pour une réaction exothermique (Q,; > o) dans un 
fluide chaud (T. > 0,) ou pour une réaction endo- 
thermique (Q}; < o) dans, un fluide froid (T. < 65), 
il est négatif dans les cas opposés, pour une réaction 
endothermique dans un fluide chaud ou pour une 
réaction exothermique dans un fluide froid. 

Quand le terme vif est négatif, il peut arriver qu'il 
compense le terme mort, toujours positif; le coeffi- 
cient de convection «. s’annule alors localement. 

Les facteurs met T.—0, dépendant de l’abscisse +, 
on ne verra pas se produire une annulation simul- 
tanée des transferts de chaleur sur toute la plaque, 
mais il pourra se faire que soit nulle la somme 
algébrique des transferts locaux de sens opposés. 

Ces conclusions éclairent certains faits qualifiés 
improprement d'anomalies, qui ont été constatés 
dans les industries chimiques sur des échangeurs 
calculés, comme à l’habitude, sans tenir compte des 
réactions intervenant dans les fluides soumis aux 
échanges de chaleur. 


20 Le profil du terme vif le long de la plaque diffère 
évidemment du profil du terme mort, puisque les 
facteurs m (27) et T. — 0, (36) qui lui sont propres 
dépendent de l’abscisse parcourue +. 

On n’en raisonnerait quantitativement que sur 
un cas concret bien défini, par tranches successives 
où l’on substituerait à l'écoulement réel une suite 
d’écoulements de Blasius tangents. 

Pour une première approximation, on peut, 
néanmoins, considérer l'écoulement réel d'ensemble. 

Dans le même esprit, nous nous limiterons pour 
l'instant au cas où l’écart maximum de tempé- 
rature T.,.—0,se maintient uniforme le long de la sur- 
face considérée, sous les influences combinées de la 
réaction et des dissociations, du rayonnement, de la 
diffusion et de convections le long d’autres surfaces. 

D'abord nul avec m au bord d'attaque (x = 0, 
x = 0), le terme vif croît et plafonne avant de 
redécroître avec lui, selon l’état de la combustion 
(facteur a), tandis que le terme mort, responsable 
par son facteur — d’un « choc thermique » au bord 

F 
d'attaque (x — o) décroît indéfiniment. 

Nous envisagerons comme exemple la flamme d’un 
gaz mixte abordant sans combustion préalable la 
surface à chauffer. 

La chaleur potentielle Q}; initialement contenue 
dans l’unité de volume du fluide pris à la tempé- 
rature initiale absolue T',, résulte pour une flamme 
de deux données pratiques fondamentales : le pou- 
voir calorifique inférieur P,; du combustible (cha- 
leur dégagée par sa combustion complète, la vapeur 
d’eau n'étant pas condensée), et le volume gazeux 
normal V;.des fumées qu'il engendre. 
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Tous deux étant rapportés à sa masse-unité, on a 
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calorifique par unité de volume normal des fumées; 
on sait qu’en combustion neutre, dans l’air, il varie 
depuis 650 gr. cal/m® pour les gaz très pauvres comme 
le gaz de haut fourneau, à 850 gr cal/mÿ pour les gaz 
riches comme le gaz de fours à coke. 


Nous admettrons . — 800 gr cal/mÿ pour un gaz 
5 


mixte. 

La densité normale p, et la chaleur spécifique 
vraie c des gaz de la flamme sont définis avec une 
approximation suffisante par 


kg/m 


et 


C= gC = (0,24 + 0,000067)£g gcal/kg masse X °C, 


T» représentant la température moyenne sur 
l'échelle considérée, c’est-à-dire ici la moyenne 
entre 8, et T.; on voit qu’il suffit de l’apprécier 
grossièrement. Elle varie depuis \800°C pour les 
fours de traitements thermiques, jusqu’à 1600° C 
pour les fours d’aciérie. 

Nous admettrons ici que 


Con= 0,32%X 1,34 g cal/m3 <0C. 
Nous aurons donc 


Q7 800 


CPo  03321,34 IRGSÈGE 


Enfin, dans les fours, l’écart thermique T, — 0, 
dépasse rarement 200°C; il peut s’abaisser à 5o°C. 
Nous admettrons ici que 


Tu — 00 = 1002 C. 


Enfin, nous poserons 


1 U. 
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SP NOEPX — = À, 
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que nous adopterons comme unité arbitraire, puisque 
nous ne cherchons que des profils en variations 
relatives. 

Il vient alors 
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Sur les figures { et 5, on a représenté, en fonc- 
tion de l’abscisse x, d’une part (courbe 1) le terme 
mort [«.l, d'autre part (courbes 2) le terme vif [al 


(fig. 4) ou le rapport .. par lequel il faut multi- 


plier le terme mort er pour obtenir le terme 
vif [ol (fig. 5). 


Les variantes 


marquées 2 correspondent à 


I * . 
=: (Ci) donc à y — x(m); les variantes mar- 


A ut \ VA 
quées 2’ correspondent à SAT donc à y — 0 


c'est-à-dire à une flamme plus « paresseuse » ou 
à un écoulement beaucoup plus rapide; les variantes 


, , « « 
marquées 2” correspondent à 10 doncày—107, 


c'est-à-dire à une flamme plus « active » ou à un 
écoulement beaucoup plus lent. 
Ces figures mettent en évidence : 


L’énorme prépondérance du terme « vif » pour 
la valeur envisagée de T, —0,;; 

Le profil à maximum de ce terme « vif »; 

La sensibilité de ce maximum (emplacement et 
amplitude) aux caractères conjugués de la flamme 
. et de l’écoulement. 


La flamme paresseuse à diffusion presque pure des 
anciens fours Martin donnera un profil tel que 2’, 
avec un maximum très retardé; la flamme brutale 
d’un chalumeau à oxygène donnera un profil tel 
que 2”, avec un maximum de [x.], plus puissant, 
très voisin de l’origine, et une valeur moyenne 
nettement plus forte, donc plus avantageuse. 

On voit que, d’une façon générale, le terme « vif » 
tend à étaler vers l’aval les effets de convection 
(courbe 2’); nous y reviendrons. L’apparence du 
choc thermique pourra néanmoins se trouver ren- 
forcée (courbe 2”), si le maximum du terme vif se 
place très près du bord d’attaque, ce qui corres- 
pond à des vitesses (u) lentes ou à des réactions (a) 

ac 


actives, PUISQUE max —= D — 1 est atteint pour 
/ max 


x-petit si - est fort. Le fait n’est pas frappant 


pour les réactions actives, qui s’achèvent prati- 
quement sur un très petit espace; il l’est pour les 
écoulements lents, par exemple dans ceux des fours 
d’aciérie récents où la flamme, lente et à combustion 
faiblement amortie, est orientée de façon à venir 
rencontrer le métal : on observe presque au point 
d'impact une fusion brutale très localisée de la 
charge, dès le début de la première phase de chauflage. 

Pour une approximation plus exacte, il faudrait 
tenir compte des variations de u, et, surtout, 


[42 : A : ’ A é A 
de - >; qui croît puis décroît pour les gaz, car a croît 
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d’abord plus vite que T, donc que u, ce qui revient 
à introduire une dilatation puis une contraction 
des tranches d’abscisses. Le maximum du terme 
vif se produit donc certainement plus tôt que ne 
l’indiquent les figures 4 et 5. 


60 


Fig. 4. — Profils de variation des deux termes du coefficient 
local «, de convection laminaire vive, en fonction du 
parcours x effectué le long d’une plaque par une flamme 
de gaz mixte offrant une température à l’infini uniforme. 


La courbe (1) représente le terme mort [z,],; les courbes (2), 


(2) et (2) représentent le terme vif [4,], pour : égal à x, 


/ 


I . A A 
TO EL UT (=). l’un et l’autre exprimés avec la même 
m 


unité arbitraire. Elles représentent également les flux calori- 
fiques [Q, l et[Q, ] exprimés avec une autre unité arbitraire. 
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Lig. 5. — Profils de variation du coefficient local a, de 


{a}, 
fe LA 


terme vif au terme mort de la convection laminaire vive, 
en fonction du parcours æ effectué le long d’une plaque 
, par une flamme de gaz mixte offrant une température à 
l'infini uniforme. 
La courbe (1) représente le terme mort [z,];, exprimé en 
unités arbitraires, ainsi que le flux calorifique [Q, h; les 


courbes (2), (2) et (2”) représentent en vraie valeur le rap- 


FA" »pour égal à 1, 1/10 et 10 (=) 
xd! 


m)°? 


COTTES- 


convection laminaire morte et du rapport local 


port de coefficients 


ainsi que le rapport === des flux calorifiques 


pondants. 
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Il faudra aussi, dans chaque cas, tenir compte 
des variations éventuelles avec x de l'écart ther- 


2000 
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Fig. 6. — Profils de variation des deux termes du coefficient 
local &, de convection laminaire vive, en fonction du 
parcours æ effectué le long d’une plaque par une flamme 
de gaz mixte offrant une température à l’infini librement 
variable avec le développement de la réaction. 


La courbe (1) représente le terme mort [a,],; les courbes (2) 
et (2) représentent le terme vif [s,], pour = () 
T,—0, = o°C et 1 000 C respectivement, exprimé avec 
la même unité arbitraire. 

La courbe (3) représente l’écart thermique total T+ — 0,en°C, 

" 


pour © — (a) T,—0, = o°C. 


mique T. — 0, qui figure au dénominateur de [x], 
Dans le cas hypothétique extrême où la tempé- 
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Fig. 7. — Profils de variation du coefficient local {æ,l, de 
[z,] 
ie 
terme vif au terme mort de la convection laminaire vive, 
en fonction du parcours x effectué le long d’une plaque 
par une flamme de gaz mixte offrant une température à 
l'infini librement variable avec le développement de la 


réaction, 


convection laminaire morte et du rapport local 


La courbe (1) représente le terme mort [«,],, exprimé en 
unités arbitraires; les courbes (2) et (2”) RE en vraie 


valeur le rapport des coefficients À pour © — I =i(E) 
1 
T,—9, = o°C et 1 0000 C respectivement. 


rature à l'infini varierait librement en fonction du 
seul degré d'avancement de la réaction, les courbes 
des figures 4 et 5 seraient à remplacer par celles des 
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figures 6 et 7, qui représentent les variations de 
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La vérité est, en général, intermédiaire et beau- 
coup plus voisine des profils figures 4 et 5 que des 
profils figures 6 et 7 


30 Le maximum du terme vif, qui est fourni par (37) 


NS I . 
où MA) atteint la valeur 


2,30 O7 
2,718 Cho( Te — 60) Ch 


[ax J2max a 


soit 
Q7 


a Le (40) 


[ax Jemax = 0,846 FT) 


Dans les conditions envisagées précédemment, 
pour un gaz pauvre, avec T,—0, uniformément 
égal à 1000 C, l’ordre de grandeur du coefficient 


‘ . : . a 
de convection vive à son maximum est pour tout =. 


1865 
To 16 [ax]4) 


[ar ,]max © 0, 846 x 
comme l'indique la figure 5. 
Ses valeurs typiques arrondies répondant aux 
limites indiquées 


È 12 , : 
650 Z V, Z 850, ‘8000 Z Ty, <1600° C 


et, 
50O0C ZT, — 0 -< 20000, 


sont les suivantes : 


Gaz pauvre, flamme très chaude, écart thermique 
maximum : 6 (x:)4. 

Gaz riche, flamme très chaude, écart thermique 
maximum : 9 (x). 

Gaz pauvre, flamme peu chaude, écart thermique 
minimum : 28 (x). 

Gaz riche, flamme peu chaude, écart thermique 
minimum : 41 (&xr)1. 


Pour les chaudières, où l’écart T,—10, atteint 
ou dépasse 10000 C, les coefficients calculés seraient 
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à diviser par 5 dans les deux premières hypothèses, 
par 20 dans les deux dernières. 

Il est entendu aussi qu’à moins de substituer 
l'oxygène à l'air, ces valeurs représentent des 
maxima qui ne sont atteints que pour une certaine 
abscisse, à laquelle («.), peut être déjà très amorti. 
Mais, dans une combustion bien réglée, les valeurs 
courantes le long de la surface ne s’en écartent que 
très progressivement quand le facteur £ =? est 
grand, c’est-à-dire quand la combustion est intense 
ou l'écoulement assez lent, comme le montre la 
figure 5. 

Multiple élevé du coefficient de convection morte 
classique, surtout pour les gaz riches et les combus- 
tions actives avec des écarts thermiques faibles, et 
soutenu sur de vastes surfaces pour les combustions 
actives dans des gaz lents, le terme de convection 
vive n’est plus, comme lui, secondaire devant le coeff- 
cient de rayonnement à haute température rapporté 
également à T.—1%, même quand l'écoulement 


gazeux reste laminaire, ce qui est assez souvent le . 


cas. 

4° L'apport du terme vif est donc en pratique impor- 
tant, parfois prépondérant, dans les fours à flammes, 
même quand ils exigent une température de régime 
très poussée (fusions, cuissons, clinkérisations, forge), 
d'autant qu’alors l’écart T',—0, est forcément réduit. 
D'ailleurs, les fusions locales observées au four d’aciérie 
ne sauraient s'expliquer par le rayonnement, qui 
n’est pas localisé, non plus que par la convection 
morte, qui n’est pas considérable. D 

Comme le terme de convection vive est lié par mau 
facteur q, ou plutôt à la moyenne de gx, on ne 
s'étonne plus d'observer quelle place la technique 
empirique des fours et des foyers fait à ce facteur q, 
sous le nom de densilé calorifique dans l’espace de 
combustion, chaque type de four étant caractérisé 
par une valeur optima de cette densité, qu’on 
exprime en grandes calories dégagées par mètre cube 
de laboratoire et par heure; par exemple : 100 000 en 
phase d’affinage et 120 000 en phase de fusion dans 
un four Martin de 6ot/h, 200 000 dans les fours 
à grille et à sole, 500 000 à 1 000 000 dans les fours 
de fusion rotatifs. On détermine très simplement le 
volume du « laboratoire » en divisant la chaleur à 
fournir par cette densité calorifique. 

Ce facteur reste cependant un peu flou, faute de 
bien préciser la fraction des calories latentes qui sont 
effectivement dégagées avant la sortie, dans des 
combustions d'activités très différentes. 

Il est néanmoins couramiment employé, surtout 
pour les projets de fours de fusion, dont on calcule 
le « cube » intérieur en divisant la chaleur horaire 
utile par le rendement cherché et par la densité calori- 
fique optima, sans qu'interviennent ni la vitesse ni 
aucun des autres facteurs des formules classiques de 
la convection, non plus d’ailleurs que ceux du rayon- 
nement ; ce qui aurait pu depuis longtemps attirer 
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l’attention des théoriciens sur l'existence d’une 
convection vive. | 

Au contraire, la considération exclusive de la 
« densité calorifique » a donné des déboires pour le 
calcul des fours à réchauffer continus, où les corps 
à traiter étant relativement froids, les écarts 
thermiques importants, les longueurs x considé- 
rables, la valeur moyenne du terme vif vis-à-vis du 
rayonnement, voire même de la convection morte, 
est fortement atténuée, ce qui rend de son intérêt 
à la considération exclusive de la vitesse. 


8. Discussion de la chaleur transmise en 
convection laminaire vive. — Envisageons main- 
tenant la chaleur échangée entre le fluide et la plaque, 
par unités de surface et de temps, qui est l'ultime 
objet des calculs calorifiques. 

Elle découle immédiatement des expressions du 
coefficient de convection; par exemple 


Qpe = CAT RE FE 80) 


, as 
— Pl oe ler. — 05) + »,30 207 | (M) 
tra Ë CPo ! 


= [Opel + [Qpele. 


Elle se compose, elle aussi, de deux termes : le 
premier représente l’apport de la convection morte, 
le deuxième l’apport supplémentaire dû à la réaction, 
que nous appellerons ferme vif. 


10 La forme du terme vif [Q,,], comporte qu’il ne 
s’annule pas avec l’écart thermique T. — 4,, ce qui 
semble paradoxal. 

En fait, il faut distinguer deux cas : 


a) D’après (36), l'écart thermique T. — 6, peut 
être nul pour toute valeur de x, parce que l’on a 


simultanément T, = 0, = o et ct 10; letralors le 
terme vif s’annule avec Q,; en même temps que lui. 


b) L’écart thermique peut être occasionnellement 
nul pour une certaine abscisse x, et une seule, parce 
qu'il y a compensation locale entre ses deux termes 
en T7, —0, et en a cette éventualité est très rare 


en pratique et exclue pour les flammes, puisqu'elle 
suppose des réactions exothermiques en fluide chaud 
ou endothermiques en fluide froid. Mais si ele se 
produit, il n’y a’ pas de raison pour que la compen- 
sation qui annule (T.—1,);, annule en même 
temps Q,,, qui fait intervenir l’abscisse x à des 
titres différents ; le fait signalé n’est donc nullement 
choquant. | 

Physiquement, il signifie que le profil T(y) à 
l’abscisse x, considérée part pour y = o de T — 0,, 
avec une tangente de pente 


(5) Ta — 00 
SN RE 0, 
DES, Ë 


puisque 7x—00=0 :(maisa#0), 
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‘pour revenir loin de la plaque à la même valeur T =6,, 
avec la même pente nulle. 

C’est donc que le profil T'(y) à l’abscisse x, présente 
au moins un maximum, et en fait, un seul. 


Fig. 8. — Profils des températures croisés au point (x, y). 


Le profil T(x) à toute ordonnée y présente lui 
aussi nécessairement un maximum, séparant les 
zones d'influence, prépondérante de la réaction et 
de l'échange simple (fig. 8). 

Cette analogie des profils T(x), et T(y)- est une 
condition nécessaire pour que soient qualitativement 
vérifiées les hypothèses que nous avons introduites, 
afin de ramener la fonction T'(x, y) des deux variables 
indépendantes x et y à une fonctionT (Z) de la va- 
riable unique Z — _ 

væ 

La part d’arbitraire que nous avons introduite 
dans la forme même de l'équation différentielle 
s’accuse (fig. 8) par le fait que T'(x), admet évidem- 
ment comme limites T,(x = o) et 0,(t — «), alors 
que-T(y). peut exceptionnellement aller de 0,(y = 0) 
à 0,(y = +). Elle se traduit dans l’équation résolue 
par l'indépendance du terme vif au regard de la 
convection morte. Désirable pour faciliter les compa- 
raisons, cette hypothèse était nécessaire pour l’inté- 
gration, et elle est admissible dans le cadre d’une 
tranche Ax. 

Sous ces réserves, le fait signalé fait comprendre 
que la chaleur créée à distance de la plaque importe 


peu vis-à-vis de la chaleur créée à son voisinage. : 


29 Le profil du terme vif (Q,.):; qui varie avec x 
comme varie 
où 
ER LR, 


CPo ÿz Vr 


part de zéro pour x = o, passe par un maximum, et 
retomberait à Zéro pour æ = «. 

Quel que soit le profil de T. — 0,, le profil de ce 
flux calorifique de convection vive [Q,,], se trouve 
être représenté à une certaine échelle par les 
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courbes (2) de la figure 4; le profil du rapport Be 


reste représenté en vraie grandeur re courbes (2) 
CAE 
Leli 
retrouve dans la courbe (2) de la figure 9. 

Si l'écart T,—06, se maintient uniforme, le 
profil du flux de convection morte [Q,,], reste 
représenté à une certaine échelle par la courbe (1) 
des figures 4 et 5. 

Si l’écart T,, — 0, évoluait librement en fonction 
de la réaction seule, le profil de [Q,,], serait repré- 
senté par les courbes (1) de la figure 9, établies 
pour une flamme de gaz mixte telle que 


. Q ren 1] 
de la figure 5, qui concernaient déjà » et qu'on 


Qr “ 
= 1865 et =: 
C20 uw 


= 2000 
+ EE eu 
w 
> Le] 
# ô 
n [ni 
5 ë 
d 1000 EL 
TO, (e] 
a . ü 
Le) 
U 2 3 < 
Fig. 9. — Profils des variations des deux termes du flux 


calorifique local Qp, de convection laminaire vive, en. 
fonction du parcours x effectué le long d’une plaque par une 
flamme de gaz mixte offrant une température à l'infini 
librement variable avec le développement de la réaction. 


Les courbes (1) et (1‘) représentent le terme mort [Qp.],, 


de x ; 
pour = égal à 1, T, —0, — o° C et 10000 C respectivement. 


a 
La courbe (2) représente le terme vif [Qp,l pour de 


et tous écarts To—,, exprimé avec la même unité 


arbitraire. 
La courbe (3) représente l’écart thermique total. T; — 6, 


a 1 - 
en °C pour I (=). D) = OC: 


Les cas courants se rapportent mieux aux figures 
ANEIREN 
Pour une meilleure approximation, il faudrait tenir 


. . . a . 
compte aussi des variations de 7’ qui pour les gaz, 


croît, puis décroît nettement. 

Provisoirement, il suffit d'imaginer des tranches 
d’abscisses dilatées puis contractées, ce qui donne 
un maximum plus hâtif. 

On s’explique donc que la flamme d’un chalumeau 
ou d’un brûleur à gaz n'offre pas son maximum 
de « pouvoir chauffant » sur la face interne du dard 
ou du cône bleu, mais un peu au delà (1 mm pour 
la flamme Bunsen), dans une zone où la combustion 
est très avancée et considérablement amortie, 
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On s'explique de même que les « tubes radiants » 
parcourus intérieurement par une flamme, empkyés 
dans certains fours modernes à « atmosphères », 
offrent un « point chaud » très marqué qui se localise 
à une certaine distance de l'émission du mélange 
combustible, distance qui est d'autant plus courte 


que 2 est plus fort. La même observation s'applique 


aux fours « poussants », qui offrent un point de 
chauffage maximum à une certaine distance du nez 
des brûleurs. L'existence de ce point chaud condi- 
tionne d’ailleurs la conception actuelle des fours 
en cause. 

On s’explique encore que, si elle offre des décro- 
chements, une voûte balayée par une flamme se 
trouve localement surchauffée après les ressauts; 
mais alors le décalage du « point chaud » est faible, 
car les remous locaux accélèrent la combustion. 


30 L'importance du terme vif dans les régions utiles 
de la flamme peut être accrue par l’amélioration 


ER 


des facteurs m et 
Co 


ac 
Êes ax a . . . 
Pour que m 7, © “ soit aussi fort que possible, 


il faut que son maximum se place dans la région 
utilisée de la flamme, donc que le contact soit 
arrêté suffisamment tôt et couvre au mieux les 
_ surfaces à chauffer; si elles sont vastes, l’écou- 
lement (u) doit être lent et la combustion (a) active, à 
amortissement lent. Ces mesures contredisent en 
partie celles qui conduisent à un bon rendement 
calorifique, comme c’est souvent le cas. Elles doivent 
être substituées à la notion trop brutale d’une densité 
calorifique moyenne. 

Pour que . = — T° soit aussi fort que possible, 
on peut agir sur le combustible, sur le comburant, et 
sur leur dosage. 

On ne gagne pas beaucoup par le choix du combus- 
tible. On démontre en effet (?) qu’en combustion 
neutre comburée par l'air, Fe varie assez peu avec 
le combustible, bien qu’individuellement P,; et V; 
varient beaucoup. Néanmoins, exprimé en grandes 
calories par mètre cube normal de fumées, il passe 
de 650 gr cal/mÿ pour un gaz de haut fourneau à 
950 gr cal/m$ pour une huile de pétrole. 

On peut donc tout de même s'attendre à ce que, 
pour une même température initiale T, des éléments 
de la combustion, l’addition de gaz riches ou de 
goudrons à un gaz pauvre augmente sensiblement 
le pouvoir chauffant de sa flamme, surtout si Ts 
n’est pas très supérieure à 0, (laquelle atteint 1550° C 
en fin d’affinage au four Martin). Ainsi s'explique 
le « doping » des fours d’aciéries, très pratiqué en 
Allemagne. 

(:) M. VÉRoN, Tendances actuelles des techniques de la 


chaleur, p.34, et Bulletin technique n° 14 de la Société Babcock 
et Wilcox. 
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On peut gagner beaucoup plus en utilisant un 
comburant plus riche en oxygène que l'air. Avec 
l'oxygène pur, V,se trouve en gros divisé par 5, 
Tr multiplié par 5, en même temps que T,— 0, se 
trouve multiplié par 2/si 0, + 6° C. 

On doit donc s’attendre à ce que le remplacement 
de l'air atmosphérique par de l'air suroxygéné 
augmente considérablement le « pouvoir chauffant » 
des flammes, grâce à l’amplification du terme vif bien 
plus qu’à l’amélioration de T. —0,, dont la possi- 
bilité reste d’ailleurs très limitée par les réfractaires 
dont on dispose. Des résultats remarquables ont été 
obtenus dans ce sens. 

Enfin, si l’on néglige les dissociations, qui n’ont 
que peu d’ampleur quand le comburant est l’air 
(T, < 18000 C), Vrest d'autant plus faible que le 
dosage du combustible et du comburant est plus 
voisin des proportions de la combustion neutre. Le 
pouvoir chauffant avec l’air est donc maximum au 
voisinage du réglage neutre. Avec l’oxygène, il est 
maximum pour un réglage nettement réducteur, 
comme le savent les soudeurs. 

On arrive aux mêmes conclusions par la considé- 
ration de la température fictive de combustion 
complète T,, qui est liée à Q. 

En effet, si l’on appelle c'e, la chaleur spécifique 
par unité de volume normal des fumées pour le 
domaine de températures compris entre la tempé- 
rature initiale T, et la température de combus- 
tion T,, le principe de l’état initial et de l’état final 
veut que 

Pai= c'en Vr(Tr— To), 


donc d’après (49) 


PTT (4) 


CPo 


T,ne dépend du combustible que dans une mesure 
assez faible; pour la combustion neutre dans l’air 
froid, l’acétylène étant excepté, elle s’étage entre 1850 
et 21500 C. 

Néanmoins, on obtient à cet égard une amélio- 
ration sensible en évitant de dégoudronner un gaz 
pauvre brut, ou en injectant du goudron de houille 
dans un gaz pauvre; son influence sur le coefficient 
de convection vive s’ajoute à son influence connue 
sur le coefficient de rayonnement, résultant du 
meilleur pouvoir émissif de la flamme. 

Mais T; dépend beaucoup du comburant. On 
l’améliore dans une large mesure en substituant à 
l'air de l’oxygène ou de l’air enrichi en oxygène. 


Enfin, T, dépend du dosage, qui doit être neutre 
avec l’air, réducteur avec l'oxygène. 


4° Considérant maintenant l’ensemble des deux 
termes de la chaleur Q,. échangée, fournie par la for- 
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mule (41), tout se passe comme si, dans la formule 
classique représentée par le terme [Q,.}, on aug- 
mentait ou diminuait l’écart thermique T7, — 0, 
réel d’un écart fictif 


mQy 


C20 


A(T a — Vo) = 2,3 


lequel représente une fraction 2,3 m < 1 de l’échauf- 
fement #L, qui (pour les liquides) résulterait d’une 


action Co plète sans échanges calorifiques. 

Ou encore, multipliant ses deux termes par cf, 
l'échange global Q,. est augmenté dans la même 
proportion que le serait la chaleur sensible c,(T;—0,) 
nécessaire pour hausser la température de 0, à T., 
si on lui ajoutait la fraction efficace mQ; de la 
chaleur latente de réaction. | 

C’est sans doute la recherche de cet avantage 
qui vient d'inspirer la conception de chaudières où 
la flamme se développe dans des tubes de fumées. 
Cette rupture avec la tradition nous semble cepen- 
dant peu indiquée, car la combustion sera certai- 
nement incomplète au contact du métal trop bien 
refroidi. 

Il faut encore voir ici la source des ennuis qu’on 
“éprouve du fait des réallumages dans les faisceaux 
de surchauffeurs, et l’un des bénéfices qu’on éprouve 
en adjoignant des brûleurs auxiliaires aux fours 
électriques à résistances. 

L'apport de chaleur supplémentaire ne s’annulant 
pas avec T. — (0, et l’emportant — sauf au début 
et à la fin du processus —., il faut s'attendre à ce 
que l’échange global Q,, ne croisse pas aussi vite que 
l'écart thermique et soit même relativement peu 
sensible à sa valeur, tant qu’elle n’est pas consi- 
dérable (supérieure à 4oo° C). Or, pour qu'il en soit 
autrement si la température 0, des corps à chauffer 
est élevée, il faudrait que la température T, de la 
flamme soit elle-même très élevée, ce qui offre 
des difficultés et des dangers, tant pour les pièces 
elles-mêmes que pour les réfractaires, et une impos- 
sibilité totale au-dessus de 18000 C dans l’état actuel 
de la technique. 

Cette vue s'accorde avec la tradition empirique 
des fours, qui limite les écarts de température T. —0, 


PUR , Do EUR ; 
généralement à es soit 150 à 2000 C au maximum; 


si la flamme est trop chaude et dépasse la valeur 
T.— 0, voulue, on la refroidit par un excès d’air. 


En confrontant les deux termes de la chaleur Q,, 
échangée, on constate que les calories apportées par le 
mélange combustible ne sont pas équivalentes quand 
elles se présentent sous la forme sensible (ces Ts) 
ou sous la forme latente (Q;), c’est-à-dire suivant 
qu'elles sont demandées au préchauffage (T,) ou au 
pouvoir calorifique (P.;). 

Les calories sensibles, n’améliorant que T. — 0,, 
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sont moins efficaces que les calories latentes, qui 
0, et TRS 
C£p0o 

Et, même, si la température 0, est forte et 
condamne l'écart T.—0, à une faible valeur, 
voire à une valeur négative vers l’entrée de la sur- 
face à chauffer, on ne peut compter sur le premier 
terme, alors que le deuxième est toujours disponible. 
A l'opposé, si la température 0, est faible et permet 
des écarts T,— 0, importants, on pourra se con- 


améliorent à la fois T, — 


. tenter d’un deuxième terme modéré. 


Un accroissement de la pression est au contraire 
également favorable aux deux termes de la chaleur 
transmise Q,, pour une vitesse u, donnée, puisque 
dans (41), les facteurs p, p, et Q, croissent comme la 


pression, qui améliore donc le facteur commun Vd 


Q7r 


sans influencer le facteur DE qui conditionne l'apport 


supplémentaire de la Sat e vive. 


9. Détermination empirique de la « cons- 
tante » de réaction a. — Pour préciser les profils 
et les valeurs moyennes de « et de Q, dans chaque 

cas concret, il faudra connaître les valeurs de a (T), 
qui dépendent du type de combustion (comburant, 
pression, degré d’homogénéité du mélange gazeux à 
l’origine). 

Si l’on a des relevés de la composition du fluide 
sur des écoulements semblables, on pourra demander 
la valeur moyenne de a sur un intervalle donné à la 
relation suivante (*?) 


\ 


qui fournit 


Puisqu'ils n’interviennent que par leur rapport, 
n et N, pourront être remplacés, s’il s’agit de 
flammes, par les teneurs relatives en gaz combus- 
tibles, qui sont faciles à mesurer. 

Pour un gaz pauvre contenant 30 pour 100 de CO 
brûlé avec son volume d'air, on aura N, — 0,15 
à 10000 C par exèmple pour l'entrée et nm, = 0,05 
à 17500C par exemple pour la sortie; avec une 
durée de parcours { — 2,5 sec, valeur normale dans 
un four de fusion, on obtiendra sur le parcours total 
qe OO de 


Du 1275710 
ETS OBS RSR PRE à 


PRE 
200 0,05 


108101, 89 © 0,3. 


De proche en proche, on pourra déterminer des 
valeurs locales de a qui permettront de localiser 
exactement le maximum du «pouvoir chauffant» Q,., 
ce qui présente un intérêt pratique évident. 


(*) On ne peut recourir à la relation (6) en q, qui donne une 
indétermination quand on la résout par rapport À a. 
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10. Extension à la convection laminaire avec 
diffusion. — Les résultats précédents sont restreints 
à un fluide où la chaleur apparaît ou disparaît 
par le fait d’une réaction, le terme complémentaire 
du coefficient de convection étant positif ou 
négatif suivant le cas. 

Nous aurions pu traiter d’une façon analogue les 
autres modes de création interne de la chaleur. 

Par exemple, dans le cas d’une diffusion spécia- 
lement intense de fluide chaud ou froid suivant la 
direction Oy dans un écoulement à deux dimensions, 
nous aurions écrit d’après (5) 


et une théorie analogue nous aurait conduits au 


coefficient de convection 
22% 
(ge) 
07 | re) 


CPO(T x — 00) 

Il est semblable à celui qui résulte de (34), K étant 
un facteur purement numérique, d’ailleurs différent 
de 2,30, car la vitesse u ne figure plus au second 
membre de l’équation (30), ce qui exclut le facteur t: 
de l’intégrale numérique J, — KJ.. 

Dans l'hypothèse mixte d’une réaction accom- 
pagnée d’une forte diffusion, cas très important en 
pratique, puisque c’est celui des flammes résultant 
du contact de deux veines gazeuses laminaires non 
mêlées au préalable, les formules précédentes (34) 
et (43) se combineraient. 

Faute de place, nous ne développerons pas pour 
l'instant ces calculs, préférant traiter plus complè- 
tement les cas de principe essentiels de l'écoulement 
laminaire pur et de l’écoulement turbulent. 


Mise 
ax = 0,332 59 VE |, Re 


IV. — La convection vive d’un fluide actif 
s’écoulant en régime turbulent permanent. 


1. Différences avec les cas précédents, et 
caractères d’une solution possible. — Dans 
l'écoulement turbulent, la chaleur est essentiel- 
lement transportée par les mouvements désordonnés 
des particules du fluide à l'échelle sensible, et non 
plus par les mouvements désordonnés des molécules 
à l’échelle atomique, comme dans l’écoulement avec 
diffusion. 

Si par exemple on envisage un écoulement à 
deux dimensions dirigé suivant l’axe Ox le long 
d’une plaque normale à l’axe Oy, le phénomène de 
la tutbulence entretient d’incessants transferts de 
particules d’une couche d’ordonnée y, dans une 
autre d’ordonnée y, < y,, et comme la température 
moyenne n’est par la même dans les diverses couches, 
on a affaire à un apport de chaleur (positif ou 


{ 


Lo 
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négatif) par transfert direct de fluide chaud ou 
froid. 

On ne connaît pas de solution mathématique définie 
pour l'écoulement hydrodynamique turbulent; il 
serait donc illusoire de former une équation de la 
chaleur, puisqu'on ne pourrait pas expliciter ses 
termes pour l'intégrer. 

Comme il est courant dans l’étude de la turbu- 
lence, on formera, par des raisonnements directs, des 
relations entre la force de frottement et la quantité 
de chaleur transférée; on envisagera les diverses 
hypothèses qui furent successivement proposées 
pour la convection morte, épousant progressivement 
la réalité déjà très complexe de ce phénomène, sans 
cependant l’atteindre exactement. 

À fortiori, la théorie ébauchée ici pour la convection 
vive ne prétend pas à la rigueur absolue. 


2. Relation fondamentale entre la chaleur 
transmise et la force de frottement. — Si le 
fluide engendre ou absorbe lui-même de la chaleur, 
due par exemple à une réaction, telle qu’une combus- 
tion, les transferts de masse d’une couche à l’autre 
auront, au point de vue calorifique, deux effets : 


1° Un transfert direct de la chaleur sensible qui 
préexiste dans le fluide en vertu de sa température : 
c'est la convection morte classique ; 


20 Un transfert de matière qui, moins usée par la 
réaction, dégagera ou absorbera dans la nouvelle 
couche plus de chaleur sensible que si elle y avait été 
prise : ce sera l’apport supplémentaire de la convec- 
tion vive. 


La notion de « parcours de mélange » imaginée 
par Prandtl permet d'exprimer facilement ces deux 
effets et de les relier au frottement. 

En chaque point (x, y) du fluide, nous désignerons 
encore par T la température, par u la vitesse projetée 
suivant Ox, et par 6 la proportion des moléeules 
initiales qui n’ont pas encore réagi; cette proportion 
décroît vers les parois, où le fluide séjourne plus 
longtemps. 

Soient (fig. 10), dans une tranche Ax, deux couches 
d’ordonnées y, et y: > y,, distantes de /, parcours 
de mélange; elles sont caractérisées par des tempé- 
ratures T,, T, et par des proportions f,, B:. 

Soient v' la composante de vitesse turbulente 
projetée suivant Oy, p la masse spécifique, c la 
chaleur spécifique de l’unité de masse. 

Par unité de surface et dans l’unité de temps, 
on admet que la turbulence transfère de la couche y: 
à la couche y, —y, — l une masse de fluide pv’, 
portant une quantité de chaleur sensible égale à 


MOT JE COT 
D a AGE (44) 


le signe étant +, car Q, est dirigée vers la paroi 
quand T,— T, est positif. 
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Soit + Q} la chaleur potentielle totale que dégage 
ou qu’absorbe la réaction complète des molécules 
initialement contenues dans l’unité de volume. 


Fig. 10.. — Schéma des notations visant les transferts 
par turbulence. 


Par unité de surface et dans l’unité de temps, on 
admettra que la turbulence transfère de la couche y; 
à la couche y, un volume de fluide v’, et ce transfert 
s'accompagne à la hauteur y, d’un apport ou d’un 
enlèvement de chaleur potentielle latente égal à 


| Q=te' Qy(B—p)=t 0 IQ = 2 (pv D SAS, (A) 


le signe étant +, car Q, est dirigée vers la paroi 
quand 8, —$, est positif. 

Sous la réserve qui sera exprimée plus loin, cette 
chaleur latente se transforme intégralement en cha- 
leur sensible créée ou perdue par unités de surface 
et de temps dans la couche y,; sinon, :le fluide non 
combiné s’y accumulerait malgré la diffusion, ou en 
disparaîtrait; $, varierait et le régime ne serait plus 
permanent. 

Le groupement (ov'l), qui entre dans les deux 
fractions (44) et (45) de la chaleur transférée sur 
l'épaisseur /, entre également dans l'expression du 
frottement par unité de surface au niveau y,, qui 
résulte, comme on le sait, d’un transport de quantité 
de mouvement : 

Soient ü, et u, les vitesses tangentielles aux 
niveaux y, Et Yo. 

Par unité de surface et dans l'unité de temps, 
on admet que la turbulence transfère à la couche y, 
une masse de fluide pv', apportant dans cette couche 
la quantité de mouvement pv'(u;, —u,) dirigée 
suivant Ox, qui tend à entraîner la couche y, dans 
le sens du mouvement. 

La résultante des forces agissant sur la couche y, 
devant être nulle en régime permanent, cette quantité 
de mouvement équilibre au niveau y, l'effort tangen- 
tiel + par unité de surface, qui est donc égal et 
orienté en sens opposé; on a 


nl 
Fe po (ua) (po Te: (46) 
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Le rapprochement des équations (44), (45) et (46) 
fournit une série d’égalités qui se combinent par 
addition 

One UT M OISTIORRENRE 


Q: pt 
(wl)= = = 
Le FN QUE C 
} 


ox , Of du 
LS 
ep dy AY.TT LP OP 0ROY 


On en tire une relation fondamentale entre la 
chaleur sensible transmise latéralement et le frot- 
tement développé tangentiellement 


LE + 0708 or QE 
CR AU HO SD 
= T1 soit encore je Ou 

dy 0y 


+ (47) 


Le premier terme figurant au numérateur du 
second membre est classique; il correspond à la 
convection morte. Le deuxième est nouveau; il 
correspond à l'apport supplémentaire de la convection 
vive. 

Dans les trois raisonnements successifs qui nous 
ont conduits à cette relation fondamentale, des 
transferts de chaleur sensible, de chaleur latente et 
de quantité de mouvement depuis la couche y, 
jusqu’à la couche y, ont été pris en considération; 
or, ils s’'accompagnent certainement de transferts de 
mêmes natures orientés en sens inverse, depuis la 
couche y, jusqu’à la couche y. 

Ces transferts en retour, qui par différence auraient 
fait apparaître des différentielles du deuxième ordre (5) 
n'ont pas été écrits, car nous admettons, comme 
on le fait implicitement toujours, qu’ils sont propor- 
tionnels aux premiers. L'erreur qui en résulte revient 
à changer la définition de {, ce qui est sans impor- 
tance puisque ce facteurs’élimine des formules finales, 
comme on l’a vu. 

La proportionnalité étant admise pour les quan- 
tités de chaleur, celle des quantités de mouvement 
ne se trouve d’ailleurs vérifiée que dans la mesure 
où le profil des températures s’identifie au profil des 
vitesses. 

Bien entendu, c’est à la chaleur latente échangée 
par ces deux transferts en sens inverses que doit 
s’appliquer la remarque faite sur la transformation 
en chaleur sensible. 


3. Expressions du coefficient de convection 
vive, pour diverses approximations faites sur 
l'écoulement hydrodynamique. — a) Adoptons 
d’abord le point de vue ancien de Reynolds-Stanton, 
qui négligent l’existence de la couche-limite lami- 
naire interposée entre la paroi et la veine turbulente. 

Évaluant grossièrement, à travers toute l’épais- 
seur de la veine, les facteurs dérivés par rapport 


() Il serait anormal que les équations différentielles 
exactes de la convection turbulente ne fussent pas du second 
ordre, comme le sont celles de la convection laminaire. 
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OT dB du  : 
y, DA D 0) qui entrent dans (47), nous obtenons 
la relation 
EE 
SN —— = ; 
CT 


Cfr Us — Uo 


LU = 9, 


Q 
C 


(48) 


donc 9, = 0. 


avec 


Q,. représente la chaleur totale transmise à la paroi 


et f- l'effort de frottement total à la paroi, rapportés 
à la surface unité; les indices o et o se réfèrent 
respectivement à la paroi et à une couche très 
éloignée. 

En fait, la deuxième fraction de Q,, représente 
une chaleur latente, mais elle se transforme inté- 
gralement en chaleur sensible sur la paroi, où le 
fluide à vitesse nulle laisse à la réaction tout le 
temps de s'achever, de sorte qu’on n’introduit pas 
d’approximation en ajoutant cette chaleur latente à 
la chaleur sensible directement transférée à la paroi. 

En se reportant à sa définition, le coefficient de 
convection totale convenant à l’approximation de 
Reynolds est 


1, 


a 
PRET ar 7e + ee | = {axli+[axle. (49) 


Te des 

Le premier terme, (x;),, représente le coefficient 
de convection morte; il doit être complété par le 
terme vif (x;),, qui représente l'apport de la réaction. 
Cette juxtaposition, déjà constatée antérieurement, 
se reproduira toujours par la suite. 


b) Adoptons maintenant l’approximation plus pré- 
cise de Prandil, qui admet l’existence d’une couche- 
limite laminaire interposée entre la paroi et la veine 
turbulente, mais la suppose à la fois mobile pour 
l’écoulement hydrodynamique, avec u(y) linéaire, et 
immobile pour l’écoulement thermique, en ce sens 
qu'il y-assimile la convection à une conductibilité 
uniforme, avec T'(y) linéaire. 

Cette correction est surtout nécessaire pour les 
liquides très visqueux, dont le nombre de Stanton 
s’écarte beaucoup de 1; elle est secondaire pour les 
gaz biatomiques et les flammes. 

La chaleur doit traverser successivement la couche- 
limite et la veine turbulente. 


1° Soient 0, l'épaisseur de la couche-limite, T; la 
température à la limite de la veine turbulente. Si 
nous admettons que la réaction y est achevée en 
raison du séjour prolongé qu'y fait le fluide, nous 
avons 


je 
Qpz = FEU me AE 


L’épaisseur à, de la couche laminaire est liée au 
frottement sur la paroi par l'équation de Stokes, 
qui fait aussi intervenir la viscosité y et la vitesse 
de transition u, 


: Ui — Ko 
fois 


AMC C/O, 
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On en déduit 


AU: 


De 


d’où la valeur du rapport caractéristique déjà pris 
en considération dans (48) ‘ 


on = et Qpx = 


Qra = S Lau (30). 
Cf UZ 
où S représente ici encore le nombre de Stanton. 
Dern 
CI 


\ 


2° Pour la veine turbulente, on aura encore, 
comme dans l'hypothèse précédente (a), 


f. a 
Dos Di) 


(Tr = TN) Ke ({ 
Œ 
— ; avec 


Q: 


es 


He OU (CE) 


Un — U 


Q, et f, comme T,;, B, et u,, se rapportent ici à la 
limite de la veine turbulente. Mais, si le mouvement 
est bien permanent et à deux dimensions et si la 
combustion s’arrête à la couche-limite, Q, s’identifie 
avec Q,. et f, avec fx. 

Le rapprochement de (50) et (51) fournit alors 
une série d’égalités qui se combinent par addition 


(Te Ti) + EL (pp) 


US Ur 


(eo £ LL (Be) 


Us + Ui(ST1— 1) 


En se reportant à sa définition, le coefficient de 
convection totale convenant à l’approximation de 
Prandtl est 


ee Be — Br Qy 
JET — 
nie Ce Joy eee ee (39) 
TES Us ts) 


La formule (52) peut être précisée de façons difté- 
rentes, suivant l'expression plus ou moins conforme 


à la réalité qu’on adoptera pour le rapport 2 ; on 


rar (44 . . 
prend en général _ — 0,5. Rien ne prouve qu'il est 
l 


ce] 


uniforme suivant Ox. 


c) Pour serrer de plus près les résultats expéri- 
mentaux relatifs à la convection des fluides dont le 
nombre de Stanton est très inférieur à 1 (huiles), 
M. G. Ribaud (#) à distingué dans la couche-limite 
laminaire deux tranches : L'une, au contact immédiat 
de la paroi, se comporte comme la couche de Prandtl; 


(‘) Nouvelle expression du coefficient de convection de la 
chaleur en régime d’écoulement turbulent (Journal de Physique, 
janvier-mars 1941, p. 12). 
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la vitesse u et la température T y croissent linéai- 
rement. L'autre, interposée entre la première couche 
et la veine turbulente, sert de transition; son profil 
de vitesse, rectiligne, prolonge celui de la première 
couche, qui, ainsi, vient se raccorder tangentiel- 
lement à celui de la veine turbulente; son profil de 
température, incurvé, se raccorde tangentiellement 
au profil linéaire de la première couche, maïs ren- 
contre le profil de la veine turbulente d’une façon 
qui échappe à la théorie actuelle, sans aucun doute 
parce qu’elle n’éprouve pas le besoin d'exprimer les 
termes difjérentiels du second ordre; le raccordement 
tangentiel ne rendant pas compte des faits, on 
opère empiriquement. 

M. Ribaud a ainsi obtenu, visant la convection 
morte, la formule semi-empirique 


Qps Cfx L 


Er = EE] —— 
k 15 < M 
| RER Fos 11) 


(33) 


qui convient bien à tous les fluides déjà expérimentés. 
Pour identifier cette formule (53) à la formule (52) 
privée de son terme vif, on voit qu'il faut poser 


TE EE 


Si donc on considérait la formule de M. Ribaud 
comme l'une des formes possibles de l’approxi- 
mation de Prandtl, elle ne pourrait se justifier théo- 


: : ; u 
riquement, car on n’imagine pas que le rapport nt 


qui ne dépend que de l’écoulement, fasse intervenir 
le nombre $ où figurent la conductibilité thermique 
et la chaleur spécifique. 

En fait, M. Ribaud a surtout voulu exprimer que 
le profil des températures diffère nettement du profil 
des vitesses quand le nombre de Stanton du fluide 
s’écarte de 1 (5). 

Mais, dans cette hypothèse, la théorie serait 
profondément modifiée et compliquée; elle devrait, 
comme nous l’avons dit, faire intervenir des termes 
difiérentiels du second ordre, ce qui explique la 
nécessité de faire appel à un raccordement empi- 
rique des profils T' (y) propres aux différentes couches. 

C’est donc en nous plaçant également à un point 
de vue empirique que nous sommes conduits à écrire 


Qpe_ 2 Ga Te ce. 


A ——— —= — 
“ T's — 0 > - 22 
IHO.5ÂS —1 


(54) 


d) Pour avoir une meilleure approximation, pre- 
nons aussi en considération l'existence de la convection 
vive dans la couche laminaire au contact de la paroi; 
elle s’apparente évidemment au cas de l’écoulement 
aminaire, traité précédemment. 


(5) Ils sont toujours un peu différents, même en régime 
laminaire, si le fluide ne réalise pas la condition S = 1, 


Reprenant l'exposé du paragraphe b, il conviendra 


: | mQr 
donc de majorer dans le rapport 12,30 CT =) 
la quantité de chaleur Q,,— Ur — 1),) qui traverse 


la couche laminaire. ARE 
Tout se passe comme si la conductibilité 4 était 
accrue suivant la même proportion, qui, dans la 
suite de notre raisonnement b, se fixe sur le nombre de 
% À 
Stanton S$ — —: 
Cu 


\ 


On posera donc 


| : 0 ou 
Su —- « RE ER Ph Aa La) 
S 1+ CETTE (33) 


4. Essai d’explicitation de la fraction restée 
disponible de la chaleur latente. — il faut 
maintenant expliciter le facteur fB, qui représente 
la fraction non encore combinée de la chaleur 
potentielle- Q;. 

Désignant par n le nombre actuel et par N, le 
nombre initial des molécules contenues dans l'unité 
de volume, on peut, comme on l’a vu, écrire d’une 
façon absolument générale 

Lt 
n To _. ARE 


Qu — 
RENTE 

En première approximation, la turbulence n’a 
pas à intervenir ici explicitement, puisqu'elle ne 
s'exerce pas à l’échelle moléculaire. Son action réside 
dans les échanges de matière entre les couches, dont 
il a été tenu compte. 

La température étant assez homogène dans la veine 
turbulente, parallèlement à la direction Oy, nous 
supposerons, pour le calcul de ce terme correctif, 
que, dans chaque tranche Ax donnée, T et a(T) 
sont uniformes, et nous désignerons maintenant 
par T une valeur moyenne dans cette tranche. 

Il vient alors 


T 
( LEA 
5 = —e at, 


a étant une « constante » caractérisant la réaction 
considérée. 
On voit que B est une fonction décroissante amortie 
du temps f. 
Si l’on admet qu'ici encore, on a 
— AS 
(12 
l’action retardatrice de la turbulence s'exprime 
par des valeurs de u plus faibles, donc de { plus 
grandes, pour x donné. 
Au total, le facteur de chaleur latente prend la 
forme 


Il décroîft donc quand x croît et quand u décroît, 


+ 
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Sous les hypothèses faites, le facteur ie COLA A 


représente la température moyenne dans une tranche, 
se mettra en facteur commun dans la diffé- 
rence B. — 6, et, dans l’expression (52) ou (54) du 
a, 


coefficient de convection, il viendra multiplier 


où p se réfère également à la température DR RARE 
Or, suivant une remarque déjà utilisée à propos du 
régime laminaire, on 4 


17 
Ces ET 


: (@ 
On peut donc remplacer le facteur variable — par 


le facteur presque constant - L. et le facteur B. —f, 
par le facteur 


SPA 
g pe ( Une um) (56) 


qui ne contient plus . ce qui apporte une double 
simplification. 


9. Expressions finales du coefficient de con- 
vection turbulente vive. — Après ce dernier 
ajustement, le coefficient de convection vive turbu- 
lente s'exprime en définitive par les formules sui- 
vantes : 


1° Si l’on adopte le point de vue de Prandt]l 
——— 
PS (Be E B )x Qr 


FE Cfx GARE Ter — 6 CPo 
PUS T+0,5(8 7 T— TE) 


= [axhi+{acl|. (87) 


2° Si l’on adopte le poïnt de vue de M. Ribaud 


tn SL EE LA Es B' Jr Qr 
re = [ox ]1 + [ox Jo à (58) 
jf 1407857 F1) 


Rappelons que :’ et B’, — f sont fournis respec- 
Qr 
P0(Z'x — Vo) 
s’explicite comme en convection nat, par (39) 
pour les flammes. 


tivement par (55) et par (56), et que = 


Rappelons d'autre part que /f. représente l'effort: 


de frottement par unité de surface de la paroi, 
connu par le coefficient de frottement 


pou (59) 
3 FU 


en fonction du nombre de Reynolds; pour un plan 


Que 
LL 


Rite = 


Sur une plaque lisse, on a, d’après Prandtl et 


Karman, 
Cp, = 0,058R,%*, (60) 


d’où l’on déduit dans ce cas particulier 


0,20 | 
fx= 0,029puiRy = 0,029 (E) pos0u80, (61) 

Notons en passant que, d’après la formule (58) 
et l'expression (59) ou (60) du coefficient de frot- 
tement, le nombre invariant 9 introduit par 
Margoulis dans l’étude de la convection calorifique 
prend la forme 


Xr fx æœ 2 
Pre GEL a ue Re7 LR= 0,029 Rue (62) 
ie Eos 
avec Rr= RE PAIE 


RE (y 
res es 1) 


6. Discussion du coefficient de convection 
turbulente vive. — On voit que, dans toutes les 
hypothèses, le coefficient de convection est formé 
ici encore d’une somme de deux termes : 

Le premier (xx), qui reste seul quand Q; =, 

est classique. Sa valeur a été calculée pour la 
première fois par les auteurs que nous avons cités. 
Nous le nommerons ferme de convection morte, par: 
opposition au second (x:), qui sera le terme de 
convection vive, pour rappeler qu’il est par exemple 
dû à une combustion vive, dans une flamme. 
. Les expressions (57) et (58) du coefficient de 
convection vive turbulente sont remarquablement 
homogènes à celles (33) ou (34) du coefficient de 
convection vive laminaire. 

Elles inspirent encore des commentaires d'intérêt 
pratique du même ordre : 

Le facteur $’ variant pour x donné dans le même 
sens que la vitesse longitudinale u, qui décroît au 
voisinage de la paroi, la différence f, —f qui entre 
dans le terme:vif est toujours différente de zéro et 
positive, en sorte que le terme vif est du signe 
de +(T.,—0,), suivant que la réaction dégage ou 
absorbe de la chaleur, comme dans la convection 
laminaire et il fallait le prévoir. 

CU 
TS CS CPo 
est plus forte ou moins forte que celle de son terme 

mn 


L'importance relative du terme vif + 


= CNT one oi a 
TRE + 0 € Po 


homologue + 
naire, suivant que B, — 82m : il est à prévoir 
qu’elle sera plus faible, la turbulence ayant évidem- 
ment pour effet d'augmenter le gradient de tempé- 
rature (5) au voisinage de la paroi, qui condi- 
fe J=0 

tionne le coefficient de convection «x. 

Mais il faut y ajouter l’influence de la combustion 
dans la couche laminaire, qui se traduit dans le 
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terme en S’ figurant au dénominateur du coefficient «.. 

Par ailleurs, $, décroissant d’abord plus vite 
que &, en raison du freinage à la paroi, puis variant 
peu, tandis que B. continue de décroître, la difté- 
rence 8. — B, qui apparaît dans le coefficient de 
convection (57) ou (58) passe par un maximum 
après un certain parcours x. L’effort de frottement 


: (BR CRE en : 
à la paroi /; = “pu décroît, lui, constamment. 


On peut donc prévoir que le terme vif (x:), du 
coefficient de convection totale doit passer par un 
maximum après un certain parcours, alors que le 
coefficient de convection morte («.). décroît toujours, 
en convection turbulente comme en convection lami- 
naire. 

Les dissemblances de la convection vive et de la 
convection morte se retrouvent donc dans tous les types 
d'écoulement des fluides, et les remarques faites au 
sujet de la convection laminaire valent encore, qualita- 
tivement, pour la convection turbulente. 

C’est ce qui nous a autorisés à émettre par antici- 
pation, à propos de la convection laminaire, des 
conclusions d’ordre technique visant des applica- 
tions qui, en fait, relèvent autant et même plus du 
régime turbulent. 

Par contre, le passage d’un régime d'écoulement à 
l’autre s’accuse nettement, en convection vive comme en 
convection morte. 

Ainsi, le frottement f. étant sensible à la rugosité 
en régime turbulent, alors qu’il ne l’est pas en régime 
laminaire, le pouvoir chauffant d’une veine fluide ne 
devient sensible à la rugosité qu’au-dessus d’une 
certaine vitesse d'écoulement, pour les flammes 
comme pour les gaz inertes, ainsi que l’a observé 
M. Meker dans ses fours à flamme hélicoïdale. 


V. — Importance technique de la convection vive. 


L'existence de la convection vive, son ordre de 
grandeur et son profil particulier soulèvent un grand 
nombre de commentaires pratiques rendant compte 
d'observations empiriques qui n’avaient pas jusqu'ici 
reçu d'explications satisfaisantes, bien qu’elles jouent 
un rôle considérable dans les diverses techniques 
de la chaleur. 

Après celles qui ont été signalées en leur lieu, citons 
les suivantes : 


1° Pour les foyers, il a été reconnu que la tempé- 
rature et la conservation des parois sont liées au 
nombre g de calories dégagées en moyenne par 
mètre cube de foyer et par heure, chaque type de 
paroi se caractérisant par la valeur maximum de gq 
qu’il peut tolérer. Cette limite est sans rapport 
direct avec les températures atteintes par la flamme, 
lesquelles ne varient pas de 10 pour 100 en valeur 
centigrade quand q varie de 300 pour 100. 


2° Sur les grilles, on a constaté que la fusion 


des barreaux métalliques se produit toujours dans 
les zones de combustion intense, favorisée par une 
distribution défectueuse de l’air (déformations de 
la grille, collage des mâchefers, soufflage de fines), 
sans que la chambre de combustion soit pour cela 
spécialement plus chaude. 


3° Dans les chaudières à écran d’eau, la vapori- 


sation croît, heureusement, plus vite avec l'allure 


de la combustion, et le rendement décroît moins 


vite, que ne le voudraient les lois du rayonnement 


et celles de la convection classique combinées. 
Cependant, l'importance de l’écart thermique entre 
flammes et parois fait qu'ici le terme vif de la 
convection est relativement faible. 


4° Dans les chaudières semi-tubulaires où les 
fumées effectuent trois parcours, d’abord autour des 
bouilleurs, puis dans le faisceau de tubes et, enfin, 
autour du réservoir, les gaz combustibles qu’elles 
contiennent éventuellement peuvent se réallumer au 
sortir du faisceau, ce qui met en jeu la sécurité de 
la plaque tubulaire bien plus que ne le commande la 
température toujours modérée des gaz en cet endroit. 


5° Dans les surchauffeurs, le moindre réallumage 
des fumées au contact des tubes chauds accentue 
la surchauffe dans des proportions deux fois plus 
considérables, malheureusement, que ne le voudrait 
la seule élévation de température des fumées. 


6° Dans les fours de recuit et de trempe traitant 
des pièces qu’il faut chauffer de façon homogène et 
modérée, dans les fours de verrerie à pots, etc. 
on doit essentiellement éviter ou retarder la 
convection des flammes, alors qu’on doit la rechercher 
et l’entretenir dans les fours de fusion des aciéries, 


“bien que dans tous les cas les gaz brûlés issus des 


flammes offrent la même température et la même 
vitesse qu'elles, avant tout échange calorifique. 


7° Dans les fours de forge continus, tels que les 
fours « poussants », où les pièces à traiter cheminent 
à l'encontre des gaz brûlés, la convection vive offre 
un grand intérêt en fin de chauffage, dans la zone où 
la température du métal, ayant atteint en surface la 
valeur convenable, reste cependant à égaliser en 
profondeur. 

Une flamme peu active, mais très chaude, compa- 
rable à celle des fours de fusion (17500 C), est bien 
capable de faire passer dans la matière traitée le 
flux calorifique qu’exigerait son chauffage à cœur 
jusqu’à la température de travail (13500 C), grâce 
à un écart thermique important; mais la tempéra- 
ture cherchée n’est atteinte en profondeur qu’au prix 
d’une surchauffe superficielle, qui engendre d’indé- 
sirables effets d’oxydation ou de fusion, notamment 
sur les arêtes. 

Une flamme très active, mais à température 
modérée (14000 C), immédiatement voisine de celle 
cherchée en profondeur (13500 C), est, elle aussi, 
capable de fournir à la matière traitée le flux calori- 
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fique qu’exige son chauffage à cœur, grâce à un 
coefficient de convection très élevé ; et la surchauffe 
superficielle est rendue impossible, puisque la flamme 
elle-même ne dépasse pas une température dange- 
reuse. Ceci implique que le profil des températures 
est plus uniforme en profondeur, donc que la chaleur 
totale entrée dans le métal est moindre. 


8 Un petit four « chauffe » toujours mieux qu'un 
gros four géométriquement semblable : c’est qu'il 
n'utilise que la zone à haute densité calorifique de 
la flamme, laissant même souvent échapper des gaz 
en combustion, ce qui compense et au delà l’aggra- 
vation des pertes par les parois, laquelle tend au 
contraire à mieux refroidir le laboratoire. 


9° Dans les récupérateurs de fours, l’admission 
accidentelle de gaz encore en combustion surélève 
anormalement la température du fluide réchaufté 
et malheureusement aussi celle des parois, qui sont 
rapidement endommagées. 


10° Une voûte épaisse ou isolée, un lit de briques 
ou un panneau radiant sont portés à des tempé- 
ratures très élevées quand une veine gazeuse en com- 
bustion active les balaie à grande vitesse, sans qu'il 
soit besoin pour cela que les réfractaires exercent 
sur la température de combustion une action cataly- 


tique, comme on l’a supposé. 


119 Dans les fours à cuve, hauts-fourneaux, cubi- 
lots, la lame gazeuse ascendante en combustion, qui 
cherche le voisinage des parois où la résistance est 
moindre, contribue à réaliser une forte capacité de 
transmission de la chaleur à ces parois. 


120 Dans les gazogènes de traction, on évite un 
échauffement exagéré des parois en écartant d'elles 
toute lame gazeuse en combustion, comme le per- 
mettent des tuyères radiales à tirage renversé plon- 
geant dans la masse du combustible jusqu’en son 
centre. 


130 Dans les flammes des chalumeaux de soudure 
ou d’oxycoupage, la substitution de l’oxygène à l'air 
joue bien plus encore en dotant la convection vive 
d’une « densité calorifique » intense par mm® x sec 
qu’en relevant de 10000€ l'écart thermique avec le 
métal. On doit diriger sur les pièces à traiter la 
pointe du dard, siège de la combustion active, et 
cependant le panache, où la combustion s'achève 
graduellement, offre assez loin une température et 
une vitesse axiale encore comparables. 


14° Dans les moteurs à déflagration, où la tranche 
gazeuse en combustion à un instant donné est un 
voile très mince si le mélange carburé est homogène 
et bien dosé, la chaleur cédée aux parois est relati- 
vement faible, malgré la haute température qu'offre 
la masse totale des gaz brâlés, leur forte pression 
et leur vitesse appréciable; elle est plus forte si la 
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combustion se poursuit derrière la flamme; elle 
devient localement intense si des «noyaux détonants » 
créent sur les parois des flammes locales sous très 
forte pression : il en résulte alors parfois des effets 
de fusion, malgré l’intermittence du phénomène (f). 


15° Dans les turbines à gaz, si la veine gazeuse 
achève de brûler au contact des aubages, elle 
leur impose des flux convectifs sévères et les met 
rapidement hors de service (6). 


16° En pyrométrie de contact, il est sensible que 
l'erreur due au rayonnement du témoin placé dans 
une flamme est inférieure aux corrections qu’on 
calcule en introduisant les valeurs connues du 
coefficient de convection morte dans l'équation 
d'équilibre calorifique du témoin. : 


17° En calorimétrie, les « calorigraphes » basés 
sur l’échauffement d’un couple par les fumées des 
gaz étudiés sont plus sensibles que ne le voudraient 
les inégalités de leurs seules températures de 
combustion, parce que les hydrocarbures, consti- 
tuants « riches », brûlent plus tardivement. 


189 Dans la propagation des incendies, une veine de 
gaz chauds offrant une température donnée (12000 C 
par exemple) et animés d’une vitesse donnée engen- 
drent évidemment des effets calorifiques plus bru- 
taux s'ils brûlent (convection vive) que s'ils sont 
éteints (convection morte). 


On citerait encore des conséquences visant les 
analyseurs de gaz combustibles à filament chauffé, 
.:2s feux de cheminée, les bacs isothermes des indus- 
tries chimiques, la balistique intérieure, la circu- 
lation sanguine, etc. 


Quand un fluide est échauffé ou refroidi, non en 
vertu de transformations chimiques, mais à cause 
de changements d’état physique, les théories précé- 
dentes s'appliquent encore, sous condition que la 
phase naissante soit très dispersée et n'’influe pas 
trop sur le frottement. 

Ainsi s’expliqueraient pour une bonne part toute 
une autre série de phénomènes, par exemple l’accrois- 
sement considérable du coefficient de convection 
qu’on observe quand un liquide devient le siège d’une 
ébullition, ou quand une vapeur devient le siège 
d'une condensation, phénomènes physiques qui 
absorbent ou dégagent une grande quantité de 
chaleur latente dans l'unité de volume. 


(5) Les phénomènes en cause ici s’accompagnent d’impor- 
tantes variations de pression, qui sont exclues des expres- 
sions actuelles des coefficients de convection vive ou morte; 
il faudra en tenir compte dans les généralisations ultérieures. 


Manuscrit reçu le 15 octobre 1942. 


GANALISATIONS A ÉTRANGLEMENTS SUCCESSIFS 


Par M. JEAN VILLEY. 


Sommaire. — Lorsque l’on fait croître progressivement le débit d’air (supposé permanent au cours 
de chaque expérience) dans une canalisation comportant deux étranglements en série, la condition 
qui localise l’apparition de la vitesse sonique dans l’un ou l’autre de ces étranglements dépend du 
rapport entre leurs deux sections, de l’augmentation d’entropie dans l’écoulement qui va de l’un à 
l’autre, et des apports de chaleur dans ce même écoulement. 

On étudie surtout ici le cas d’un écoulement adiabatique. La vitesse sonique ne peut apparaître dans 
l’étranglement d’amont que si sa section s, est plus petite que celle s, de l’étranglement d’aval, avec 
un écart relatif d'autant plus accentué que les décoordinations entre s, et s, sont plus importantes. 
Pour que l’étranglement d’aval puisse devenir lui aussi sonique, il faut que le rapport « 


Sa 


ne dépasse 


pas une certaine limite d’autant plus grande que le rapport — Se de la section maxima entre les deux 


étranglements à celle du premier étranglement est plus grand. Si le rapport : 


So 


$ est plus grand que cette 


2 


limite, l’onde de choc passera brusquement, au cours de l’augmentation progressive du débit, de la 
position S, vers une position en aval de l’étranglement s,, avec variation discontinue de la vitesse 
d'écoulement en s4, qui passe d’une valeur subsonique à une valeur supersonique. 

Le paragraphe 9 ajoute quelques indications relatives aux modifications que des échanges de 
chaleur localisé entre les deux étranglements apportent à ces résultats. 


1. Écoulements adiabatiques. — Lorsque, dans 
une canalisation parcourue par un courant gazeux 
et qui présente plusieurs étranglements successifs 
en série, on augmente progressivement la diffé- 
rence des pressions entre les deux extrémités pour 
faire apparaître des vitesses soniques locales, la 
succession des phénomènes observés variera sui- 
vant les dimensions relatives des divers étrangle- 
ments. 

Pour éclairer la question, nous étudierons en 
particulier dans quelles conditions la vitesse d’écou- 
lement peut être simultanément égale à la célérité 
du son dans deux cols situés en série sur un même 
écoulement d’air permanent. 

Nous envisagerons l'hypothèse de l’adiabaticité 
rigoureuse de l’écoulement, et nous nous conten- 
terons ensuite d'indiquer succinctement au para- 
graphe 9 les modifications qu'introduisent des 
échanges de chaleur avec l'extérieur. 

Dans le cas de l’adiabaticité (et sans qu'il y ait 
lieu de supposer l’irréversibilité, qui donnerait 
alors l’isentropie), la température T d’un col sonique 
est déterminée par la température T, de l’état initial 
de départ. La: vitesse u en un tel col satisfait en 
effet simultanément aux équations : 

HET OR UIE (vitesse sonique), (1) 


et 
= C(To— T') 


(conservation de l'énergie), (2) 


qui donnent ensemble 


YRT =2C(To5— T), 


d’où 
RE To , of 2 
TEE — $ 
‘Re TO AT ONE 
CEADOURE Nr 4 
Te US Su (97 


? 


Alors la vitesse d'écoulement est la même aux 
deux cols soniques dont nous représenterons les 
sections par 8, et 54. L'égalité des deux débits (équa- 
tion de continuité) exige donc que les densités soient 
inversement proportionnelles aux sections, autrement 
dit que les volumes spécifiques v, et v, soient propor- 
tionnels aux sections 5, et 5. Mais, sur l’isotherme T, 
les volumes spécifiques et les entropies varient 
dans le même sens. Or, le long de l’écoulement 
adiabatique qui va de s, à 54, il est impossible 
que l’entropie diminue; elle peut augmenter s’il 
y a des irréversibilités (frottements ou ondes 
de choc). La conclusion, c’est que le col sonique 
d’aval peut avoir une section s, plus grande que 
le col sonique d’amont s,. L'hypothèse inverse 
est impossible. , 

On peut donc affirmer que si, dans une canali- 
sation, supposée adiabatique, qui comporte plusieurs 
étranglements en série, celui dont la section est 
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minima est situé en aval des autres, seul il pourra 
devenir sonique. , 

Pour préciser notre étude, envisageons le cas 
simple où la canalisation est constituée par deux 
tuyères de Laval disposées en série à la suite l’une 
de l’autre. 

Appelons S; l'entrée du premier convergent 
et & son col, S, la section maxima dans la région 
qui est simultanément la fin du premier divergent 
et l’entrée du second convergent, s, le second col, 
et enfin S, la section maxima du second divergent. 


FIST 


Supposons que S, constitue la sortie d’un réservoir 
d’amont À et que S,; débouche dans un réservoir 
d’aval B. 


2. Écoulement par succion en aval. — Consi- 
dérons d’abord le cas où l’état P,, T, est maintenu 
invariable dans le réservoir d’amont; par exemple 
en le laissant en communication très large avec 
l’atmosphère. 

On peut alors provoquer et entretenir l’écou- 
lement en réalisant par pompage une dépression 
dans le réservoir d’aval B. Nous appellerons P,4 
la pression qui y est entretenue. 

La température T, n’a d’ailleurs pas d'influence 
sur l'écoulement. Cette particularité est liée à la 
très faible vitesse de la conduction thermique, 
qui ne peut pas remonter l’écoulement; alors Te 
ne peut pas agir sur celui-ci, tandis que la propa- 
gation des perturbations de pression avec la vitesse 
du son leur permet de remonter l'écoulement tant 
que sa vitesse reste subsonique : la pression P, 
réagit donc en amont. 

Mais, dès qu'il s'établit, en s, ou en s4, un col 
sonique suivi d’une petite portion, si courte soit-elle, 
d'écoulement supersonique, la pression P,, qui 
peut encore agir sur la partie de l'écoulement en 
aval de la région supersonique, n’a plus aucune 
action en amont de celle-ci : ses variations ne modi- 
fient donc plus le débit. 

Des considérations simples d’homogénité (1) per- 


G) Cf. C. R. Acad. Sc., 214, (février 1942), p. 210. 
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mettent de conclure que la loi de variation du débit 
ést de la forme 
Po 


UT, (4) 


en appelant k le rapport k — _ 
0 


D = f(x) 


de la pression 


d’aval à la pression d’amont. 

Il est clair que, lorsque l’on fait décroître k progres- 
sivement, à partir de la valeur 1, la fonction f (k) 
va en croissant. Elle devient ensuite constante 
lorsque X devient inférieur à une certaine valeur h 
pour laquelle l’écoulement devient sonique dans 
une des sections étranglées s. La valeur limite À 
est déterminée, pour un état donné d’amont (P;, To), 
par la valeur de la section sonique s, si l’écoulement 
est isentropique jusque là. 

Le débit sonique est en effet alors : 

1 
DE? Po 


an 2 di STE . 
D=s(—) Ve 7e: #2 


et l’on a la condition 
1 


po=s( 2) 


qui détermine la limite h, et qui permet de la 
calculer car on connaît la fonction f (k) dans l’hypo- 
thèse où l'écoulement reste pratiquement isentropique 
jusqu’en S;. On a alors en effet 


Y—1 


Ts te P3\ * P3\028571 # 
lis (A () 


et la condition 
ui = 2C(To— Ta) = >CTa(— 7.) 


0 
donne, en notant que, dans l’écoulement subso- 


nique en S;, on a P; = P,, et en posant : 6 
0 


dy V AE co (8) 


d’où, pour le débit, l'expression : 
Ji 
D Ps Saus = Sses Vert y ). 
Mais l’isentropie donne aussi 
1 

Ps Y & 
D — KŸ— 0,71#28. 
Pc ; (9) 


d’où il résulte finalement : 


D = S3po VE) 


Le Po 
ou, en utilisant p, — AT; 
Ps V2C ( lo 
or 2 ee ES (10) 
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La fonction f (k) est donc : 


FR) =Ss Vers nl A 


. (11) 


pour 1 > k >> h; et la limite h est déterminée par 
la condition (6) qui donne : 


1 


Ro) a ant) 
£ T— ï = 12 
Se R ( £ (—) R(y +1) (52) 


ou, 


| DNS ne 2YR 
542857 511428 0,5 — | —_ ——— 
MERS EEE 


R 1} pt » 
ou enfin, en observant que Te Dee donne à l’ex- 


. . 1 
pression sous le radical la valeur & 


s 
(A1:42857 2 P1,71428 0,5 — ee es 


S- (3) 


—@— = 0,2588 


Lorsque X devient inférieur à cette valeur À, la 
fonction f (k) devient constante avec la valeur 


1 


SN enr te cUus 
k =s ) Vre Ne ss 
J(X) Y+1 ROSE UR 651,2 7 


Elle donne dans (4) pour expression du débit sonique 


Es = $.0,68477 ae 
VPo%o P0 


La succession des phénomènes lorsque k continue 
à décroître au-dessous de À va différer suivant 
que le col sonique s’est amorcé en $, ou en sg. 


(14) 


0,5 
Des 
— —— 


(15) 


3. Col minimum en aval. — Lorsque s, est 
plus petit que s,, la remarque générale faite au 
paragraphe 1 nous indique que le col sonique qui 
s’amorce en s, restera unique. Le débit restant 
constant malgré l’abaissement progressif de P,, 
l’écoulement ne deviendra jamais sonique en &. 

La diminution de P, ne modifie que la partie 
de l’écoulement située en aval de 5,, et il est facile 
de voir la nature de la modification progressive 
qu’elle lui impose. 

Tant que l'écoulement reste subsonique dans 
la section de sortie S;, la condition P, — P, reste 
satisfaite. Donc P, diminue progressivement, alors 
que la pression P,; au col sonique reste invariable- 
ment égal à 0,527 P, (si l’on admet que l’écoule- 
ment est pratiquement isentropique depuis À 
jusqu’en $4). 

Cette diminution de P, est liée au déplacement 
de l’onde de choc qui sépare la partie supersonique 
de l'écoulement, en aval immédiat de 54, de la 
partie restée subsonique, en amont immédiat de S. 
L'augmentation globale de pression (P,; — P,) com- 
porte : une diminution de P, à P, dans la partie 
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supersonique du divergent, une augmentation brusque 
(P,—P;) dans l’onde de choc, et une augmen- 
tation de P, à P, dans la partie subsonique du 
divergent. À mesure que la section S; qui contient 
l’onde de choc se déplace de sç à S,4, cette somme 
algébrique, qui mesure (P4— P;) diminue comme 
on peut le vérifier par des calculs faciles, mais 
assez longs (?). 

Pour une certaine valeur h' de k, l’onde de choc S, 
arrive dans la section terminale S,. À partir de là, 
c’est-à-dire pour k <h', la pression P;3 dans la 
section de sortie cessera d’être égale à la pression P, 
dans le réservoir. Elle subit une brusque diminution 
au moment où l’onde de choc traverse la section S3, 
et se fixe à une valeur invariable p, que l’on calcule 
facilement si l’on admet que l'écoulement est 
isentropique jusqu’en Ss. 

L’équation de continuité 


pe Su = const., (16) 
et l'équation de l’enthalpie 
= 2C(To TT) (AT) 


jointes aux conditions d’isentropie, écrites en posant 


12 : Ne 2e 
— = % pour alléger l'écriture, soit : 


Po 
# 
o) — £0 rT— 20 7x 0,71428 (18) 
et 
Y—1 
= Tofn \ "TT 0:28574% (18°) 
donnent : 
poT°:71428 S /2 CT (1 — 70:28571)0,5 — const. 
ou 
S(r112856 — m1,71427 )0,5 — const. (19) 
ou 
S2(71,12856 —. rl11#27) — const. (192) 
Maison a 
Lai) 
2 Y—1 fre NS 
T6 = ( ) == (=) 10020: (20) 
Y+I 1m) 


On admet en général plutôt la valeur 0,527. 
En réalité, on ne saurait prétendre avec certitude 


à cette précision, car le rapport y =£ n'est pas 


constant lorsque la température varie, et comporte 
des incertitudes de l’ordre de 3 millièmes : la valeur 
que l’on admet comme la plus exacte dans les 
conditions courantes est 1,405 plutôt que 1,4 que 
nous avons gardée parce que beaucoup plus simple. 

Des évaluations très simples des dérivées loga- 
rithmiques montrent que l'erreur relative sur 
est nettement plus grande (et de sens opposé) 
que l'erreur relative sur y, et que les erreurs rela- 


(2) GC. R: Acad. Sc., 1942, 214, pp. 345 et 872. 
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tives sur les exposants de x dans l'équation (19) 
sont du même ordre de grandeur (et de sens opposé) 
que celle sur *. 

On ne peut donc pas compter sur les trois derniers 
chiffres inscrits dans les exposants de (19), et 
nous Calculerons simplement 


128 


1, 4 
0,527 0,527 = 0,067065. 


Gardant la notation x pour représenter r4 — Le; 
l’équation (19') devient * 


S8(rit28— x1,714) — 5,2 0,067065, (21) 
qui détermine r et permet de le calculer par substi- 
tutions et interpolations. | 

Lorsque S; est beaucoup plus grand que s,, 
cette équation a une racine très voisine de : et 
une autre très voisine de zéro. La première corres- 
pond à l'écoulement isentropique partout subso- 
nique (), la seconde à l'écoulement isentropique 
dont la partie s$— S; est supersonique c’est 
cette dernière racine qui nous intéresse ici et que 
nous désignerons par t. 

Pour k < h' il ne se produit plus aucune modifi- 
cation dans la canalisation. 

Tant que k > r, on passe de p, à P, en traversant 
des ondes obliques de choc situées dans la veine 
libre. En réalité, par inertie, la pression dépasse P, 
et subit des variations oscillatoires réalisées, avec 
variations de section de la veine, en traversant 
alternativement des ondes obliques de choc et des 
ondes obliques de dépression. 

Lorsque k devient plus petit que 7, c’est-à-dire 
lorsque P, devient inférieur à p4, on aura d’abord 
dans la veine libre des ondes obliques de dépression, 
suivies de variations oscillatoires comme dans le 
cas précédent. 

Mais ces phénomènes complexes, localisés dans 
le réservoir B ne retentissent pas dans la canali- 
sation, où l’écoulement est subsonique de S;, à s, 
(avec un maximum de vitesse en $,, mais sans y 
atteindre la vitesse sonique) et supersonique de 5 
à S3g. 


4. Cols égaux. — Si les deux cols $&, et s; sont 
égaux, la vitesse y deviendrait sonique simulta- 
nément dans l'hypothèse de l’isentropie parfaite. 
En réalité, par suite des décoordinations par frot- 
tement, c’est inévitablement le col aval s qui 
deviendra sonique. 

On a en effet 


P;vÏ > Pav} ou FT en ER € À) 
qui exprime l’augmentation d’entropie, et 
s S V6 po c 
us= us Qu M —— (23) 
6 pe LA Us 


(3) Avec un maximum de vitesse en ss. 


# 
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qui exprime la conservation de la masse; d’où 
Tuto = Tour, (24) 
Mais 
ui? = 2C(To— T5) et u»2= 2C(To— To) 
permettent d’écrire (24) sous la forme 
Le des 
To To— T5) * > Ta(To— To) *? . (24) 


Si l’on fait croître le débit de telle sorte que T, 
et T; s’approchent simultanément de la valeur 
sonique 


2, ml To 
EP ee 
on écrira 
T : 
Te = ko d et TESTÉ ee 
m2 A0) 


où k, et k, sont tous les deux supérieurs à 1. Il 
s’agit de déterminer lequel atteint le premier la 


valeur 1. 
La condition (24’) devient : 


(24) 


L'hypothèse k, — 1 donne alors 
0,202 > k:(1,2 — ko)? 


ou, en posant k = 1 +e, 
(1+e)(0,2—e)?< 0,202 


ou 


0,2 
(1 + eo, (1— =) «40,202; 
\ 


0.2 
que l’on peut écrire, lorsque e est très petit 


? € 
Ge) (io 2) <i ou I—E<I. 
50 


Cette condition est toujours satisfaite, c’est-à-dire 
er , D 

que T}, peut devenir égale à us alors que € est encore 

positif, donc que T, est resté plus grand que cette 

valeur sonique _ 


2 

Au contraire l'hypothèse k, — 1 conduirait à 
la condition k, (1,2 — k)** > 0,2°** qui est impos- 
sible puisqu'elle aboutirait à 1 — € > 1. 

C’est donc inévitablement le col aval s, qui devient 
sonique, et, à partir de ce moment, si l’on continue 
à faire baisser la pression P,, on provoquera, en 
aval de 54, les modifications successives décrites 
au paragraphe ©, sans qu’il en résulte aucune modi- 
fication au col & qui reste subsonique avec un débit 
constant. 

Pour que la vitesse devienne sonique simulta- 
nément en & et s il faut donc en réalité une sec- 
tion s; un peu plus grande que la section s,, soit 
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Se = S(1 +») avec n positif et d'autant plus 
grand que l’augmentation d’entropie entre & et 54 
est plus importante. 

Mais, même en supposant que l’on atteigne 
simultanément T, — Lie CLR _ donc la 


vitesse sonique aux deux eridroits, c’est le col 5 
seul qui évoluera pour devenir ce que l’on appelle 
normalement un col sonique, c’est-à-dire un col 
où la vitesse, subsonique en amont, devient super- 
sonique en aval. En effet, s’il se formait en aval 
de $, une région supersonique, et corrélativement 
une onde de choc, l’augmentation d’entropie de & 
à s4 irait en croissant à mesure que cette onde de 
choc s’éloignerait de s,. La température T, restant 
la même, le volume spécifique v, irait en croissant, 
et cela ne serait possible que si l’on faisait croître 
la section s, proportionnellement à w. 

On a donc là, une fois réalisé le bon rapport (1 + ») 
entre les sections s, et 5, le moyen de maintenir 
de façon tout à fait stable (4), malgré les variations 
éventuelles de la pression P, du réservoir aspirateur, 
une vitesse qui passe en $&, par un maximum égal 
à la célérité du son. Si l’on avait en s, une section 
plus petite que s, on pourrait y réaliser aussi une 
vitesse passant par un maximum sonique, mais à 
condition de maintenir dans le réservoir aspirateur 
une pression P, rigoureusement invariable : si la 
dépression aspiratrice diminue, la vitesse en 5 
devient subsonique; si la dépression aspiratrice 
augmente, la vitesse u, reste sonique, mais elle 
n'est plus maxima et l’on a des vitesses super- 
sonique, entre s; et l’onde de choc qui s'établit 
un peu en aval. 


5. Col sonique en amont. — Supposons que & 
soit nettement plus petit que s,, de telle manière 
que la vitesse y devienne sonique alors qu’elle 
est encore subsonique en 5,. 

Si l’on continue à faire baisser P,, la pertur- 
bation mécanique se fait sentir jusqu’au voisi- 
nage de s,, à travers s, où la vitesse passe par un 
maximum, mais reste subsonique. Cette pertur- 
bation est une diminution de la pression, et elle 
provoque, à l'aval immédiat de s,, la région super- 
sonique que l’on avait prévue, au paragraphe 3, 
en aval de 54. | 

De cette diminution de la pression P, nous ne 
pouvons plus tirer de conclusion immédiate relative à 
la variation de la température T,, et corrélativement 
de la vitesse u;, car nous n’avons plus la relation 
d’isentropie entre s, et sç (puisque l’on traverse l’onde 
de choc S;); mais l’invariabilité du débit, main- 


. REX u 
tenue par le col sonique & entraîne que 7 reste 
constant .dans une section donnée quelconque, 
et par conséquent en particulier en aval de S', 


(#) Si P, est invariable. 
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et nous pourrons déduire le sens de variation de u 
du sens de variation de v. Or l'augmentation progres- 
sive d’entropie produite par cette onde de choc dont 
la discontinuité augmente, entraîne que PvY va en 
croissant; cela donne 


PyeY-1de + #Y dP > 0. 


La pression P diminue, en même temps que P:, 
dans tout l'écoulement subsonique en aval de S',, 
c’est-à-dire que dP <o; d’où résulte dv > 0. Le 
volume spécifique v va donc en croissant, et par 
conséquent aussi la vitesse d'écoulement u. Dans 
la section minima 5, où la vitesse est maxima, 
u, pourra atteindre la valeur sonique si l’augmen- 
tation d’entropie est suffisante, c’est-à-dire si le 


> d i se produira si °* 
rapport -* est assez grand, ce qui se produira si 


S3 _, : è 
est plus grand que la valeur : : nécessaire (celle-ci 


2 


va en croissant en même temps que la section 5). 


6. Deux cols soniques simultanés. — Suppo- 
sons ainsi atteint le régime sonique en 5. On a 
alors 


Us —= Uo et 


L'égalité des débits (5) donne ; 


On peut alors calculer l'augmentation (S;—S,) = (S)5 
de l’entropie en utilisant l'expression 


S = cLogT + R Loge + const., 


qui donne simplement ici, puisque Ty = T2, 


S6 


(S)S = R Log © = R Log © 

Si l’on suppose négligeable l’augmentation d'en- 

tropie par frottements dans les écoulements en 

amont et en aval de l’onde de choc, il faudra iden- 

tifier cette expression (25) avec la variation 

(S', —'S3) — (S);' réalisée par la discontinuité de 
l’onde de choc. On a 


7 0 


(SX = clos + R Log (26) 
et l’on obtient l’équation 
R Loë SP LE A RDS (27) 
D PAT 1817 =) DSe “ 


Si l'on suppose s, donné et si l’on cherche la 
position S; où doit alors être l’onde de choc, celle-ci 


(5) Le débit en *, n’est plus donné par la relation (5), 
parce que l’écoulement n’est pas isentropique jusque-là, 
du fait de l’onde de choc; c’est 5, qui règle le débit par la 
formule (5). 
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s’introduira sous la forme de sept inconnues S;; 
Us, Vs Ta; Us, V,, T, liées les unes aux autres par 
six autres équations, et l’on se trouvera en présence 
d’un système très complexe de sept équations à 
sept inconnues simultanées. 

Le problème devient infiniment plus simple si 


on le prend en sens inverse, c’est-à-dire si l’on 


cherche à calculer, par l’équation (27), la section 
sonique $ç qui correspond à une position S; donnée 
de la première onde de choc, car on sait alors évaluer 
_ les diverses caractéristiques us, vs, Ts; u', 0, T\. 

En effet, u, v, et T, sont connus, et, quand on 
se donne S:, l’état 3 à l'entrée de l’onde de choc 
est déterminé par les équations : 


SE 2 9 Uo nl ù 
_ —= nl (continuité), (28) 
3 D > 
Toi Top (isentropie), (29) 
uÿ —uÿ = 2C(To— Ts) (enthalpie). (30) 


D'autre part l’état 3’ à la sortie de l’onde de choc 
est calculable, à partir de l’état 3, en écrivant les 
équations de continuité, de l’enthalpie et de la 
quantité de mouvement, qui donnent : 


ur = CET, TS), (32) 
Be Gu—u)= PP R(TE TT). (39) 
P3 P3 P3 


On peut donc construire la courbe qui donne 
les valeurs s,; correspondant aux positions S, succes- 
sives de l’onde de choc. 

Si la section s; du second col correspond sur 
cette courbe à une section S; plus petite que S$S, 
le col $, devient sonique au moment où l’onde de 
choc arrive dans cette section S;, où elle s’immo- 
bilise; s, réalise dorénavant une barrière contre 
toute réaction vers l’amont des variations de P,. 
Si l’on continue à faire baisser P,, il se produit 
en aval de s, une seconde onde de choc S, qui s'éloigne 
progressivement comme on l’a vu au paragraphe 3, 
sans modification de l’onde S;. 


7. Col d’aval supersonique. — Si l’onde de 
choc parvient jusqu’à la section maxima S, sans 
que le col sç soit devenu sonique, on arrive à un 
régime instable. Il n’est pas possible que l'onde 
de choc orthogonale continue à se déplacer progres- 
sivement vers l’aval, de S, vers 5. En effet, cela 
supposerait une vitesse décroissante dans la partie 
supersonique du convergent, donc une Giscontinuité 
décroissante dans l’onde de choc, c’est-à-dire une 
augmentation d’entropie de plus en plus faible. 
Il en résulterait une augmentation plus considérable 
de la pression depuis la valeur P, (déterminée par P, 
puisque l’écoulement est pratiquement isentropique 
jusqu’en S,) jusqu’à la valeur réalisée dans une 
section donnée d’aval, par exemple P,. Cela exigerait 
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donc une pression P, croissante, alors que nous 
étudions des succions croissantes par abaissement 
progressif de P,. 

il n’y a donc qu’une solution possible, c’est que 
l'onde de choc, au moment où elle dépasse la section 
maxima $S4 se trouvera brusquement balayée 
au delà du col s, pour aller chercher une position 
convenable S, — si elle existe — dans le divergent 
final 54 Sa. 

Alors l’étranglement 5, ne réalise plus rien d’autre 
qu’un minimum de section dans l’écoulement super- 
sonique S$,S,. Il pourra toutefois rendre l’hypo- 
thèse d’isentropie nettement incorrecte dans la 
portion S,s,; de cet écoulement. Il semble en effet 
se produire en général des dissipations notables 
dans les convergents supersoniques, où se propagent 
des perturbations obliques dont la superposition 
peut donner des discontinuités. S’il en résulte une 
augmentation d’entropie non négligeable auprès de 
celle que devrait donner l’onde de choc S$, si elle 
intervenait seule pour maintenir P, — P,, on ne 
saura plus calculer la section S, nécessaire; on sait 
seulement que sa valeur initiale est alors inférieure 
à S,. Si S, est plus petit que la section S, nécessaire, 
on passera directement de l’onde de choc ortho- 
gonale en S$S, aux ondes obliques en aval de S;. 
Si au contraire l’onde orthogonale retrouve une 
position S, convenable, elle continuera à se déplacer 
vers S, lorsque P, continue à baisser, puis sortira 
en aval de S, et sera remplacée par des ondes obliques 
dans la veine libre. 

Là encore, faute de savoir évaluer l’augmentation 
d’entropie dans le convergent supersonique $,s4, 
on ne saura pas Calculer la valeur p, à-laquelle 
se fixe la pression dans S, après la sortie de l’onde 
de choc. 


8. Écoulement par compression en amont. — 
Dans le cas où'la pression d’aval P, est maintenue 
invariable — par exemple par une large commu- 
nication du réservoir B avec l’atmosphère — et 
où l’on provoque l’écoulernent en élevant progres- 
sivement la pression P, dans le réservoir d’amont A, 
les choses semblent à première vue plus complexes 
du fait que le débit ne reste plus constant malgré 
la présence d’un col sonique. En effet, même si 
l’on maintient T, invariable, ce qui maintient 
invariable la vitesse sonique, l’augmentation pro- 
gressive de P, augmente partout la densité, et par 
conséquent augmente le débit dans le col sonique. 

Cela ne modifie pas toutefois l'équation géné- 
rale (4) du débit, ni l’expression (11) de la fonc- 
tion f (k) valable, dans l'hypothèse de l’isentropie, 
tant que 7 = _ reste plus petit que la valeur " 

9 
déterminée par l’équation (13). 
Cela ne modifie pas non plus l'expression (14) 


de la valeur constante que prend f(k) lorsque _ 
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devient plus grand que 7° ni l'expression (15) du 


débit qui en résulte; la seule différence c’est que ce 
débit sonique (15) n’est plus constant, par suite 
des “’ariations que l’on impose à P,. En se reportant 
à l’équation (4) on voit que, si T, est maintenue 
constante, le débit sonique D varie proportion- 
nellement à P,. 

De même, le résultat que le col d’amont &, ne 
peut pas devenir sonique à moins d’être plus petit 
que le col d’aval s,, subsiste. 

S1 Sg < S2 (1 + n), où n est déterminé par l’augmen- 
tation d’entropie que provoquent les frottements 
entre $& et sç au moment où le régime sonique va 
s'établir, c’est en s, que s’amorce le col sonique. 

Laissant de côté les complications qu'intro- 
duisent les décoordinations par frottements, nous 
supposerons ceux-ci négligeables pour une étude 
de première approximation (n négligeable). 

Nous envisagerons donc d’abord l’hypothèse 5, <$ 
pour laquelle le col sonique s’amorce en 5. On a 
alors P, = 0,527 P,,-qui va en augmentant propor- 
tionnellement à P,. Pour que P, reste égale à P,, 
il faut que l’augmentation de pression (P3 — Ps) 
aille en diminuant : cela sera obtenü, comme au 
paragraphe 3, par le déplacement progressif d’une 
onde de choc S; dont la discontinuité augmente à 
mesure qu’elle se déplace de s, à S,. Lorsque cette 
onde vient traverser la section terminale S;, la 
pression en S;-tombe brusquement de P,=P, 
à ps < P%ÿ et l’on aura alors. des ondes de choc 
obliques dans la veine libre. Mais ici, dans la suite 
de l’expérience (où P, est maintenue constante), 
ps Va croître proportionnellement à P,. Lorsque p; 
devient égale à P,, les variations oscillatoires de 
pression dars la veine libre cessent. Lorsque ps 
devient supérieure à P,, elles recommencent, mais 
en débutant par des ondes obliques de chute de 
pression au lieu des ondes obliques de choc du 
Cas Ps < Py 

Au cours de toutes ces modifications, l'écoulement 
reste subsonique en amont de s,, et la pression P, 
dans la section $, > s, reste plus grande que P,, 
c'est-à-dire plus grande que 0,527 P,; c’est-à-dire 
que le col s, ne devient jamais sonique : si T', reste 
constant, l’augmentation du débit y est réalisée 
par augmentation de la densité sans augmentation 
de la vitesse d'écoulement. 

Considérons maintenant le cas s, < 54, avec un 
écart assez net pour que les augmentations d’entropie 
par frottements ne risquent pas de mettre en défaut 
les prévisions que nous faisons en les considérant 
comme négligeables. 

Le col sonique s’amorce alors en &,, et, si l’on 
continue à faire croître P;,, le raisonnement du cas 
précédent, basé sur l’augmentation de P, et l’inva- 
riance de P;,, montre qu’une onde de choc S$, doit 
naître en s%, et s’avancer progressivement vers 
l’aval pour introduire une augmentation d’entropie 
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croissante. Les calculs du paragraphe 6 restent 
applicables pour déterminer si, avant qu'elle ait 
atteint la section maxima S,, le col s, devient lui 
aussi sonique. 

Envisageons cette hypothèse. Nous ne pouvons 
plus alors, comme dans le cas de la succion aval, 
invoquer l'absence de réactions mécaniques sur 
la région comprise entre S, et s; pour conclure à 
une position invariable de l’onde de choc S; au cours 
de la suite de l’expérience (augmentation progressive 
de la pression P,). Mais, si nous examinons le système 
d'équations (27 à 33), nous observons que l’état 
d'amont P,T, n'y intervient pas directement, 
mais seulement par l'intermédiaire de us, To, V. 
Or uw et T, restent invariables dans l'hypothèse 
que nous envisageons où l’on maintient T, inva- 


riable, et v, n'intervient que par le rapport à qui 


est déterminé par : puisque l’on admet que l’écou- 


lement est isentropique entre 5 et S,. Par conséquent 
la position $S,; de l’onde de choc reste invariable 
dans l'hypothèse où T, est invariable et où les frot- 
tements sont négligeables entre & et S3. 

Si la section étranglée s, est trop grande, elle ne 
deviendra pas sonique, et l’onde de choc S;, au 
moment où elle arrivera dans la section maxima $,, 
sera balayée pour aller se reformer dans une section S, 
du second divergent, ou même elle sera balayée 
totalement si la section finale S, n’est pas supérieure 
à la section S, nécessaire pour rétablir P;, = Po. 


9. Écoulement non adiabatique. — Lorsque 
l’on abandonne l'hypothèse de l’adiabaticité, l’équa- 
tion de base (2) est modifiée. Si nous appelons Q; 
la quantité totale de chaleur X4 dg reçue par l’unité 
de masse dans tout son parcours en amont de la 
section S; (où les termes successifs dg peuvent 
être positifs ou négatifs, et sont exprimés en unités 
mécaniques; ils sont supposés invariables sur chaque 
élément du parcours, condition pour que l’écoulement 
soit permanent), cette équation (2) devient, en utili- 
sant P;, T;, v, et u; pour désigner les caractéris- 
tiques de l’écoulement non adiabatique dans la sec- 
tion sonique 5; 


(34) 


qui, avec 
LME SERA (35) 
donne 
YRTi=2C(T— Ti) +20, 
d’où 
SC OQ; 
RE xt 
" yR+2C (ro | 
ou 
ne Qi 
Ti = = (r+ &) (36) 


\ rs 


di à 


No 8. à 


Le débit sonique est modifié par suite des modifi- 
cations simultanées de la vitesse et de la densité. 
Sa valeur est multipliée .par 


us Pi Ti 
ui Ps Ti 
Mais, d’après (35) 


er AE 
OT HAUTE s 


d’où, d’après (36) : 


À T'; ès To 0,5 
NT CEE ER 
À és m.) 


qui est déterminé par la valeur globale Q; de la 
chaleur reçue par l’unité de masse. 
Par contre l’altération de la pression, c’est-à-dire 


(37) 


le rapport ee dépend de la façon dont les apports 
i 


élémentaires de chaleur ôg ont été distribués au 
cours de la variation de température, car la variation 
de densité dépend, pour une variation de tempé- 
rature donnée, de la variation d’entropie 


où les Ôq, outre leur intervention directe dans le 
premier terme, peuvent-avoir une action dans le 
second sur les quantités d'énergie décoordonnée, 
par exemple en déplaçant une onde de choc. 

Si nous prenons pour l’entropie l’expression 


S — CLogT — RLogP, 


nous avons alors 


S, — S;= C Log 4 — R Log. (38) 
On a d’autre part : 
LES 
L Le Ôæ 
S—S= [ po 
} 
s ds. 7? 
Se f : FE à PA 
d’où . 
dæ 
ss= f % - +[ (FF). 
PHARE ES | 


En égalant les expressions (38) et (39), et en 
utilisant (36), on obtient alors : 
ï.)" (40) 


ff (re 


Au total, compte tenue de (37), les apports de 
chaleur modifient le débit sonique dans le rapport 


( CT 1 P: 
CT + Où P; 


CTi+ Q; 
CT 


R Log = = CLog ——— 


(4) 
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P; 
P; 

Lorsque toutes les décoordinations dw de l’écou- 
lement adiabatique sont considérées comme négli- 
geables, le débit sonique est donné par l’expréssion (5) 
du débit sonique isentropique, multipliée par le 
facteur (41) lié aux apports de chaleur. On doit 
alors, dans (40), faire dw — o, et, sauf cas excep- 
tionnels (6), on aura aussi dw' — o au même ordre 
d’approximation. 

On pourra admettre cette hypothèse à titre de 
première approximation lorsqu'il n’y a pas d’onde 
de choc en amont de la section sonique, c’est-à-dire 
en particulier dans le cas du col sonique unique. 
Dans le cas des deux cols soniques simultanés, 
le débit est réglé par le col d’amont, pour lequel 
l’approximation reste valable dans les mêmes con- 
ditions. 

Alors, si nous posons 


où est la valeur déterminée par l’équation (40). 


= Her (42) 
et 
© Log 4 — | ae D. (43) 
on aura par (40) : 
. —\eDi (44) 
Le débit sonique est alors, en partant de (15) 
D!= s,0,68427 Va Je. (45) 


La condition pour que la vitesse atteigne la 
valeur sonique simultanément dans les deux cols & 
et s, n’est plus l'égalité 5 — 54 (?). Elle s’écrira 


eB: eB S° À 9 0,5 
So  — $; ou — — eBi—B; ( = ; 16 
HAS AUS S6 na HN 
ou x 
muse [ar A 0,5 
pe = A 4; 2 RT/ ( À) , (46) 
S6 6 
avec 
ACTE QE ES Cr Gr 
A» CTo+ Où CT CTo+ Q2° 
où est mis en évidence l’apport de chaleur (Q; — Q:) 


localisé entre s, et 54. 
La condition pour que le col sonique s’amorce 


en & est donc 
59 < eBi—Bs EE Fer 
ñ A 


Si cette condition est satisfaite, 


(47) 
il reste à déter- 


(£) Où un apport de chaleur provoquerait de notables 
perturbations mécaniques. 

() Dans la théorie de première approximation qui suppose 
les frottements négligeables; en réalité $ç = 8, (1 + n). 
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iiner celle pour que s, puisse devenir lui aussi sonique. 
On l’obtiendra encore en étudiant la correspon- 
dance entre la section S,; de la première onde de 
choc et la section s, du second col sonique. 
Les équations du paragraphe 6 sont modifiées 
comme suit. On a d’après (36) 


Ts OT EU ES (OR 


Te t4 CTo+ (E 


Us (Z)"= = Se)" 
& CTo+ Qù ; 


et 


d’où 


Po Sa Uo 


06 us SG Us RL CTo+ (O7 0,5 SG 
CTo+ Qù S9 


On obtient donc, à la place de (25) l'équation 


Se B= CL ge 
+ Ro,5 Log ER A +R Log® 
ou 
Re (C0 0) Log Ep : +RLog®. (48) 


L’équation (26), relative à l’onde de choc, n’est 
pas modifiée, et, en égalant les expressions (26) 
et (48) on obtient, à la place de (27) l'équation 


SE + R Log © 


(CR Res - 


T; C2 
= c Log — + R Log =. 
T; p3 


(49) 
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Pour la détermination de l’état 3 à partir de la 
donnée S:, on a : | 


— l'équation (28) de continuité, qui n’est pas 
modifiée ; 


— à la place de (29), celle qui exprime la variation 
de l’entropie par les échanges thermiques, soit 


T : ‘à 
cLog gi +RLoge = k 7 (50) 
— et, à la place de (30), l'équation de l’enthalpie 
modifiée par l’apport de chaleur 


(51) 


Quant aux trois équations (31), (32) et (33) 
relatives à l’onde de choc, dont l'épaisseur extré- 
mement faible ne laisse pas place à des échanges 
de chaleur, elles ne sont pas modifiées. 

Les quaritités de chaleur (Q: — Q) et Q,, et la 


U3? — U3? = 2C(T2— T3) + 2(Q:3— Q2). 


* dq em : ; 
somme 7 qu interviennent dans ces équa- 


tions, dépendent de la position S; considérée pour 
l’onde de choc. 

Lorsque s, est plus grand que la valeur qui cor- 
respond à S; — S,, en arrivant en $, l’onde de choc 
sera balayée en aval de $, si l’on continue à faire 
baisser P, ou à élever P,;. 


Manuscrit reçu le 15 mars 1943. 


